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1.  INTRODUCTION 


In  recent  years  high-energy  propellants  have  been  developed  containing 
cyclotetramethylenetranitramine  (HMX)  to  improve  the  performance  of  rocket 
motors.  In  developing  such  propellants,  several  accidental  explosions  have 
occurred  during  rocket  motor  firings.  The  purpose  of  this  study  is  to  iden¬ 
tify  the  causes  of  the  explosions  so  that  such  occurrences  may  be  prevented 
in  the  future.  The  problem  selected  for  study  was  that  of  pressure  build¬ 
ups  within  burning  propellant  cracks.  This  selection  was  based  upon  I  IT 
Research  Institute's  (IITRI)  experiments  with  burning  secondary  nigh  explo¬ 
sives  in  closed  bombs. 1,2  In  these  experiments  Composition  B  and  PBX  9404 
undergo  multiple  cracking  prior  to  a  violent  explosion  or  high-order  deton¬ 
ation.  Two  facts  support  the  presence  of  multiple  burning  cracks.  The  first 
is  random  firing  of  ionization  probes  placed  within  4x4  inch  cylinders  of 
the  two  secondary  explosives.  Random  firing  commenced  when  the  pressure 
reached  about  70  bars.  The  second  is  fragments  of  Composition  B  (See  Figure 
1)  found  in  the  bomb  following  sudden  venting  of  the  bomb.  Explosive 
fragments  varied  in  dimension  from  about  1/10  Inch  to  1  Inch.  All  fragment 
surfaces  were  covered  with  a  "frozen"  melt  layer  indicating  that  they  were 
burning  prior  to  sudden  pressure  relief. 

The  subject  matter  reported  herein  is  directed  at  the  three  project 
objectives.  In  order  of  execution,  these  are: 

(1)  Use  of  analysis  to  identify  mechanisms  and 
propellant  properties  causing  pronounced  pres¬ 
sure  transients  within  burning  propellant 
cracks. 

(2)  Compare  predicted  pressure  transients  with 
shock  wave  pressures  and  durations  known  to 
initiate  propellants. 

(3)  Conduct  experiments  with  which  to  complement 
and  validate  analytical  predictions. 

DDT  was  studied  analytically  in  two  phases.  In  the  first  phase,  cracks 
were  considered  to  propogate  into  a  cavity  of  high-pressure  high  temperature 
gas  such  as  the  combustion  chamber  of  a  rocket  motor.  This  problem  is 
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Figure  1.  Composition  B  fragments  found  in  closed  bomb 
folloinng  sudden  pressure  relief. 
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Illustrated  by  Figure  2a.  It  was  chosen  In  that  sudden  exposure  of  cracks  to 
high  pressure  gases  creates  stress  waves  that  subsequently  act  to  compress 
cracks  after  being  reflected  from  the  motor  case.  In  this  study  two-demensional 
models  were  developed  with  which  to  predict  gas  pressures,  velocities  and 
temperatures  within  the  crack,  anc.  the  propellant  response  to  the  predicted 
gas  conditions.  Mechanical  responses  of  cracks  to  reflected  stress,  waves  and 
transient  pressures  of  gases  within  cracks  were  included.  This  model  allows 
one  to  examine  the  consequences  of  various  cavity  gas  pressures,  cavity  tem¬ 
peratures,  motor  cases,  crack  dimensions  and  propellant  properties  upon 
pressure  build-ups  within  cracks.  This  model  provides  efficient  predictions 
of  dynamic  pressure  build-ups  to  a  few  thousand  atmospheres  or  bars.  Above 
such  pressures,  variations  of  crack  widths  and  burning  rates  along  the  length 
of  crack  became  too  great  to  afford  efficient  computations.  The  models 
greatest  value  is  In  establishing  crack/propellant  conditions  caused  by  sudden 
exposure  of  cracks  to  high-pressure  high-temperature  gases  that  may  subsequently 
promote  extremely  rapid  and  pronounced  pressure  rises  needed  to  cause  detona¬ 
tion.  Among  these  conditions  are  substantial  melt  formation  over  the  burning 
propellant  crack  surfaces,  minimal  crack  widths,  and  the  presence  of  high- 
amp!  Itude  stress  waves. 

The  second  analytical  phase  focused  upon  those  portions  of  the  crack  with 
conditions  most  likely  to  lead  to  rapid  pronouned  pressure  transients.  An 
illustration  of  this  problem  is  provided  by  Figure  2b.  Gas,  propellant  and 
crash  conditions  are  considered  uniform  alona  the  length  of  the  crack  element. 
The  latter  implies  negligible  gas  escape  from  the  crack  element.  Neglect  of 
gas  loss  is  predicated  upon  occurrence  of  extremely  transient  pressures  within 
a  few  tenths  of  a  millisecond  or  less.  Durinq  such  periods  aas  disturbances 
travel  less  than  10  cm.  This  model  was  computerized  and  used  to  determine 
values  of  the  stress  wave  amplitudes,  amounts  of  melt,  and  crack  widths  needed 
to  generate  pressures  of  the  order  to  tens  of  K  bars.  The  two-dimensional 
model  serves  to  indicate  the  validity  of  the  above  values  used  in  the  one¬ 
dimensional  model. 

Predicted  pressures  transients  varied  widely  with  the  propellant/crack/ 
stress  wave  values  used  in  the  one-dimensional  model.  These  studies  identi¬ 
fied  conditions  needed  to  generating  pressures  of  the  order  to  10  K  bars 
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Figure  2a.  Crack  configuration  treated  by  PROS  two-dimensional  model. 
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Figure  2b.  Crack  configuration  treated  by  PROS  one -dimensional  model. 
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within  a  few  tenths  of  a  microsecond.  Each  of  the  above  conditions  involved 
applyinq  high  amplitude  stress  waves  to  relatively  narrow  crack  elements  con¬ 
taining  substantial  melt. 

Validation  of  model  DDT  predictions  was  attempted  experimentally.  The 
experimental  design  simulates  partial  closure  of  a  burning  propellant  crack 
caused  by  a  stress  wave.  Unfortunately,  DDT  was  not  achieved  due  to  inade¬ 
quate  of  the  void  space  over  or  between  the  burning  propellant.  Inadequate 
closure  was  attributed  to  more  rapid  acceleration  of  the  burning  than  was 
expected.  It  is  suspected  that  the  rapid  acceleration  was  caused  by  ejection 
of  propellant  melt  or  foams  into  the  hot  combustion  gases  by  the  rapid  qas 
flow. 

The  remainder  of  this  report  contains  analytical  basis  and  predictions 
of  the  models,  and  conclusion  drawn  from  the  study  insofar  as  DDT  is  con¬ 
cerned. 
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2.  ANALYTICAL  PROCEDURES 


This  section  shall  review  means  used  to  predict 

•  gas  flows,  gas  pressures  and  gas  temperatures 
produced  within  cracks  suddenly  exposed  to 
high-pressure  high-temperature  cavities  of  gas 

•  stress  wave  amplitudes,  deformation  of  cracks  by 
the  stress  waves  qas  pressures 

•  heating  and  burning  of  propellant  crack  surfaces 

Nomenclature  is  presented  in  Appendix  F  along  with  values  used  for  the  various 
parameters. 

2.1  PROPELLANT  HEATING  AND  BURNING 

Propellant  heatinq  will,  of  course,  vary  depending  upon  whether  or  not 
the  propellant  is  burning.  Prior  to  ignition,  the  propellant  will  be  heated 
connectively  by  the  hot-gas  stream  and  to  a  lesser  extent  by  decomposition 
of  heated  propellant.  The  latter  is  termed  internal  heating.  Following  ig¬ 
nition,  the  Internal  heating  increases  and  an  additional  heat  source  comes 
into  be1ng--namely  conduction  from  the  flames.  For  the  above  reasons  each  of 
the  above  periods  will  be  discussed  separately. 

This  section  is  primarily  concerned  with  dynamic  burning  of  propellants 
such  as  HMX  with  a  melt  or  foam  layer  over  the  burning  surface.  The  presence 
of  a  foam  layer  makes  the  analysis  more  complex  than  that  of  propellants  with¬ 
out  such  a  layer.  After  reviewing  the  equations  used  for  propellants  with  a 
foam  layer,  we  shall  Indicate  revisions  needed  for  propellants  that  burn  with¬ 
out  a  foam  or  melt  layer. 

Figure  3  represents  symbols  used  to  represent  various  temperatures  (T^, 
Tm,  T),  heat  fluxes  (qf,  qp,  q)  and  regression  rates  (r^r)  associated  with 
the  melt  or  foam  layer  and  the  two  Interfaces  of  the  foam  layer. 
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Temperatures 
Heat  Fluxes 
Regression  Rates 


Figure  3.  Principal  variables  used  to  predict  propellant  burning. 


2.1.1  Heat  Fluxes 


Pn'or  to  the  start  of  melting,  the  convenience  heat-flux  q  equals  that 
entering  the  solid  propellant.  It  is  described  by  the  expression  used  by 
Kuo.’  In  terms  of  the  gas  temperature  Tg  and  the  temperature  of  the  pro¬ 
pellants  surface,  the  convective  heat  fli*x  q  is 

f  s  hc  (Tg  •  V  <» 

while  the  heat  transfer  coefficient  hc  is  given  by 


During  melting,  there  are  two  sources  by  which  the  melt  or  foam  layer  is 
heated.  The  most  obvious  is  by  convection  produced  by  the  hot  gas  stream. 
This  heat  flux  is  described  by  Equation  1.  Heat  is  also  generated  within  the 
melt  due  to  thermal  decomposition.  The  latter  heat  flux  is  described  by 

Q  •  Mf  •  Z  »  exp  ( -E/Tf )  (3) 
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where  represents  the  mass  of  melt  per  unit  area  of  surface.  The  expres¬ 
sion  Z  exp  ( -E/T^.)  Indicates  the  fractional  rate  of  decomposition  of  the  mass 
The  constant  of  proportionality  Q  represents  the  quantity  of  heat  transfered 
to  the  inelt  by  the  evolution  of  a  unit  mass  if  gas.  Lacking  experimental 
data  we  shall  assume  Q  equals  the  Internal  heat  Q$  transferred  during  burning. 

During  melting,  the  sum  qf  of  the  fluxes  entering  the  melt  is 
qf  =  hc  (Tg  -  Tf)  +  Qs  *  Mf  .  Z  exp  -  (E/Tf)  (4) 

Part  of  the  above  flux  acts  to  heat  the  foam,  while  the  remainder  is  expended 
In  melting  the  propellant  and  in  heating  the  underlying  solid  propellant. 

The  portion  that  acts  to  melt  and  heat  the  solid  propellant  is  represented 
by  qp.  In  turn,  the  portion  of  qp  that  is  conducted  into  the  solid  propel¬ 
lant  is  represented  by  q.  The  difference  between  the  heat  fluxes  qp  and  q 
is  expended  in  melting  the  propellant  so 

qp  -  q  =  pr  Qm  (5) 

where  r  represents  the  rate  of  melting  of  a  propellant  having  a  density  p  and 

a  heat  of  fusion  Q.  Expressions  for  q  and  q  are  described  below, 
m  p 

The  heat  flux  q  is  described  in  terms  of  a  heat-transfer  coefficient  h 
P 

by 

qp  =  h(Tf  -  TJ  (6) 

Clearly  h  depends  upong  foam  motions  produced  by  escaping  gases.  For 
this  reason  h  Is  considered  to  depend  upon  the  rate  of  gasification  in  the 
following  fashion: 


h  =  ci  [exp(-  i  )]C2  (7) 

'f 

where  ci  and  c2  are  constants  which  are  evaluated  in  Appendix  B  for  HMX. 
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Substitution  of  the  above  expression  for  h  into  Equation  6  yeilds  the 
following  expression  for  q^. 


qp  =  c,  Cexp(-  ^)]Ci(Tf  -  Tm)  (8) 

At  the  solid  propellant  side  of  the  melt  interface,  two  conditions  must  be 
satisfied.  The  first  expresses  the  rate  of  heat  conduction  q  Into  the  solid 
propellant  as  follows: 


q  =  -K|I  (9) 

The  second  expresses  the  constant  temperature  as  the  melt  Interface  as  follows 


T 


(10) 


In  view  of  the  dependence  of  q  upon  temporal  and  spacial  variations  of  the 
temperature  within  the  solid  propellant,  the  heat  flux  q  is  computed  numer¬ 
ically.  The  procedure  for  computing  the  temperature  T(x)  of  solid  propellant 
is  described  in  Appendix  A  for  time-dependent  fluxes  q.  The  above  procedure 
does  not  depend  upon  the  presence  of  foam.  However,  foam  will  affect  the 
fluxes  q.  When  propellants  melt  q  may  be  obtained  from  Equations  5  and  8; 
when  propellants  do  not  melt,  q  equals  q^..  The  above  statements  also  apply 
after  burning  commences.  In  this  regard,  a  discussion  of  propellant  ignition 
is  presented  in  Appendix  C. 


W zaX  Fluxea  \jtzn  StaM  oj  Burning 

Once  burning  starts,  the  heat  flux  entering  the  propellant  is 


"f  '  T7  tV VV  *  VTf  •  Tm>  *  W  +  "f  %  * 


-Bpr, 


»rf  «w  <VV  *  \  <w  t-p^r>  <»> 
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Equation  11  includes  a  cross  flow  p^u.  When  the  cross  flow  Is  negligibly 
|  small,  Equation  11  reduces  to  that  developed  by  Krler,1*  namely 

»  . 

prJ 

j  vf 'VVV'.'W'W* 

j  >rA  *  (VTf>  W) 

|  Initial  terms  of  Equations  11  and  12  account  for  heat  feedback  from  the  flame; 

the  second  term  pr^Q$  accounts  for  internal  heating  near  the  propellant  sur¬ 
face;  the  third  term  corrects  the  previous  flux  when  the  specific  heat  C  of 

i  m 

the  melt  differs  from  that  Cq  of  the  gas;  and  the  last  term  of  Equation  11 

j  accounts  for  heat  fluxes  produced  by  gas  flowing  laterally  over  the  burning 

:  propellant  surfaces.® 

i  Bars  over  variables  of  Equations  11  and  12  Indicate  their  values  during 

steady  burning  at  the  pressure  present.  Heat  fluxes  depend  upon  the  pressure 
(  p,  burning  rate  rf  and  cross  flow  pqy.  The  importance  of  the  three  sources 

I  of  fluxes  changes  as  the  burning  becomes  more  dynamic.  Pressure  rises  will 

accentuate  heat  feedback  from  the  flame  and  erosive  heating.  More  rapid 
)  burning  diminished  the  latter  fluxes  and  increases  internal  heating.  As  the 

burning  becomes  increasingly  dynamic.  Internal  heating  becomes  more  and  more 
j  dominant.  In  this  regard  over  90  percent  of  the  heating  can  be  due  to  inter¬ 

nal  heating  during  periods  of  highly  dynamic  burning  needed  to  cause  DDT. 

i 

|  In  view  of  the  importance  of  internal  heating,  it  Is  wise  to  briefly 

examine  assumptions  implicite  In  the  expression  pr^Qs  of  Equations  11  and  12, 

)  and  their  conseouences  upon  predicted  burn  rates.  First  of  all,  this  expres¬ 

sion  indicates  a  linear  dependence  between  the  flux  and  the  rate  of  gas  evolu¬ 
tion  pr^..  Such  a  relationship  presumes  that  the  amount  of  heat  transferred 
to  the  propellant  is  proportional  to  the  amount  of  evolved  gas  regardless  of 
,  its  temperature  or  its  residence  time  In  the  foam.  In  this  regard,  rapid 

i  pressure  rises  will  act  to  compress  gases  within  the  foam,  and  thereby  In¬ 

creases  the  period  during  which  the  gases  undergo  partial  decomposition  and 
j  transfer  heat  to  the  melt.  For  the  above  reason,  It  Is  believed  that  the 
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expression  Pr^Qg  used  for  internal  heating  under  estimates  the  Internal  heat¬ 
ing  during  periods  of  rapid  pressure  buildup. 

Another  presumption  of  the  expression  pr^.Q$  Is  that  the  foam  is  not 
blown  off  into  the  flames  by  the  rapid  evolution  of  gas.  This  subject  as 
well  as  the  latter  will  be  returned  to  in  Section  6. 

2.1.2  Steady  and  Noiuteady  Boon  Ratea 

Steady  burning  rates  rf  of  Equations  11  and  12  are  given  as  a  function 
of  pressure  P  by  the  conventional  equation 

rf  =  aPn  (13) 

where  a  and  n  are  constants  that  are  determined  experimentally.  Nonsteady 
burning  rates  r^.  are  a  function  of  the  foam  mass  M^.  and  temperature  as 
expressed  by  the  following  Arrhenius  relationship 

rf  =  MfZ  exp  (-E/Tf)/p  (14) 

where  Z,  E  and  p  represent  the  frequency  factor,  activation  energy  and  den¬ 
sity  of  the  propellant,  respectively.  The  expression  Z  exp(-E/Tf)  represents 
the  rate  of  gasification  of  the  foam  on  a  fractional  basis  at  the  temperature 
Tf.  An  expression  for  rf  is  developed  in  Appendix  D  for  propellants  without 
a  foam  layer. 

2.1.3  TKom  Heat,  HaM  and  TempyutfuM. 

Here  we  shall  represent  the  sensible  heat  within  a  unit  surface  area  of 
foam  by  Qf  and  reference  it  to  the  propellant  melt  temperature  T  .  Rates  of 
change  of  are  given  by 

V  Of  -  «p  *  rfC*  <VV  (15) 

Here  q^-qp  represents  the  net  heat  flux  to  the  foam  while  the  last  term  rep¬ 
resents  the  rate  of  heat  carried  away  by  the  escaping  gases. 
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■Tates  of  change  of  the  foam  mass  are  described  by 


Mf  =  p(r-rf)  (16) 

When  the  sensible  heat  of  the  foam  is  referenced  to  its  melt  temperature  T  , 

5f '  WW-  ^  »7> 

Tf“CT7*Tm  (18> 

m  f 


2.2  GAS  DYNAMICS 

As  noted  earlier,  two  computerized  models  were  developed  to  study  the 
behavior  and  effects  of  gases  within  propellant  cracks.  The  first  model  is 
two-dimensional.  It  predicts  the  temperature,  pressures,  and  velocities  of 
the  gases,  crack  deformations,  heatinq  and  burning  of  the  propellant  surfaces, 
and  the  generation  of  stress  waves.  This  model  predicts  special  variations 
of  the  burning,  crack  widths  and  pressures  alonq  the  length  of  cracks. 

The  second  model  focuses  upon  those  portions  of  the  crack  with  which 
appreciable  pressures  rises  may  occur.  This  model  assumes  uniform  conditions 
along  the  length  of  the  crack  element  and  hence  is  one-dimensional. 

The  one-dimensional  model  serves  to  indicate  crack/stress  waves  conditions 
needed  for  DDT;  the  two-dimensional  model  serves  to  identify  conditions  pro- 
noting  DDT. 

2.2.1  Two-VZmniionaZ  ttodll  oJ\  Ga&U 

Equations  used  to  describe  the  flow  and  behavior  of  gases  within  cracks 
are  presented  in  Aopendix  E.  They  are  written  with  respect  to  an  Eulerian 
frame  of  reference  as  functions  of  time  and  distance  along  the  length  of  the 
crack.  The  equations  conservation  of  mass,  momentum,  and  energy  as  well  as 
the  effects  of  inertia,  wall  friction,  crack  deformations,  and  the  addition 
of  mass  and  energy  from  the  burning  propellant  crack  surfaces. 


Rates  of  mass  Mq  and  energy  addition  from  the  burning  propellant  are: 

1^  =  2fprf,  and  (19) 

\  -  \  C«r  +  **  +  V*  +  Cm  <VV]  '  <J°> 

The  above  equations  apply  to  a  unit  length  of  crack.  The  factor  f  of  the 
above  equations  represents  the  ratio  of  the  area  of  the  irregular  crack  sur¬ 
faces  to  that  of  a  planar  surface,  while  the  factor  2  accounts  for  the  two 
crack  surfaces. 

For  steady  burning  Equation  20  reduces  as  it  should  to: 

\  ■  2fprfQr  (21) 

Reactions  of  the  evolved  gases  are  assumed  to  be  instantaneous. 

2.2.2  0ne.-Vin\en6ionaJt  ModeZ  oj  Gaj>ej> 

The  one-dimensional  model  assumes  uniform  gas/propellant  conditions  along 
a  given  length  of  burning  crack.  Rates  of  increase  of  the  internal  energy 
are  described  by 


61  -  (Q  -PC  -  M  )/(p  C  ) 

'  yy  rt  '  '  K  n  u i ' 


at 


g  w 


(22) 


where  and  are  given  by  equations  19  and  20,  respectively. 

Combustion  gases  within  the  crack  are  assumed  to  obey  the  Nobel -Abel 
equation  of  state 


P(l/pg-t»  *  RT9 


(23) 


where  the  internal  energy  l  is 

*  =  P(l/Pg-b)/(Y-l)  (24) 

where  the  ratio  of  the  specific  heats  of  the  gas  is  assumed  constant. 
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Time-dependent  pressure  rises  of  the  gas  are  described  In  terms  of  the 
Interval  energy  by 


P  =  (Y-m/(l/Pg-b) 


2.3  STRESS  WAVE  GENERATION  AND  EFFECTS 

Stress  waves  are  generated  by  the  sudden  exposure  of  cracks  to  gas  cav¬ 
ities  with  higher  pressures  than  initially  exist  within  the  cracks.  High- 
amplitude  stress  waves  can  also  be  generated  by  extremely  rapid/pronounced 
pressure  buildups  due  to  burning.  Analyses  describing  the  generation,  wave 
properties,  reflection  of  stress  waves  from  motor  cases,  and  their  effects 
upon  cracks  are  described  in  Appendix  E. 

In  the  remainder  of  this  section  we  shall  present  equations  describing 

the  consequences  of  stress  waves  of  amplitude  AP  upon  the  crack  width  C  . 

w 

The  velocities  of  the  two  crack  walls  will  be  distinguished  by  Wj  and  W2, 

where  Wj  equals  the  velocity  of  the  crack  wall  upon  which  the  stress  wave  Is 

incident.  Velocities  Wj  and  W2  are  considered  positive  when  they  act  to  in- 

* 

crease  the  crack  width  C  and  negative  otherwise.  The  velocity  Wi  of  the 

w 

crack  wall  subjected  to  the  stress  wave  is  given  by 


Wi  =  r. 


P-P  2AP 
0  - 


The  velocity  W2  of  the  other  crack  wall  is 


W »  =  r 


f  +  WT 


Notice  that  Wj  may  be  either  positive  or  negative  while  W2  Is  always  positive. 
This,  of  course,  is  because  the  gas  pressure  and  burning  act  to  expand  the 
crack,  while  the  stress  wave  acts  to  contract  the  crack. 

Rates  of  change  of  the  crack  width  Cw  are  given  by 


-  W,  +  W2 
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Initially  a  stress  wave  contracts  a  crack;  therefore  the  buildup  pressure 
p  within  the  crack  causes  the  crack  to  expand.  Gas  pressure  continues  to 
buildup  during  crack  expansion  before  decaying. 

2.4  COMPUTATIONAL  PROCEDURES 

Three  different  predictive  schemes  are  used  to  assess  the  effect  of  stress 
waves  upon  burning  propellant  cracks.  These  involve  dynamic  predictions  of 

(1)  gas  behavior 

(2)  crack  deformations 

(3)  burning  propellant 

Stable  and  accurate  predictions  require  careful  choice  of  time  steps  for  each 
of  the  above  predictions.  In  this  regard,  order  of  magnitude  smaller  time 
steps  are  usually  needed  for  Items  1  and  2  than  needed  by  Item  3.  The  latter 
is  particularly  true  during  the  early  stages  of  pressure  buildup  wherein  the 
propellants  response  to  the  changing  gas  conditions  is  relatively  slow.  For 
the  above  reasons  the  modeling  of  Items  1  and  2  was  treated  separately  from 
Item  3. 

2.5  COMPUTATIONAL  PROCEDURES 

Predictions  are  made  in  a  stepwise  fashion  with  respect  to  time.  Each 
model  uses  output  supplied  by  the  other.  Burn  predictions  Supply  rates  of 
energy  and  mass  flows  into  the  crack  for  the  predictions;  the  gas/crack  pre¬ 
diction  supply  gas  pressures  for  the  predictions. 

2.5.7  BuAn  Pxaczdwizi 

Burn  predictions  require  calculations  of  the  following  time-dependent 
terms: 

•  foam  temperatures  T^. 

•  melt  rates  r 

•  foam  mass 

•  sensible  heat  in  the  foam 
following  each  time  step 
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Time  steps  are  calculated  as  follows 


minimum 


4c*2  °-lprf 

fC 

r  :  f 


(29) 


where  r-  represents  the  mean  melt  velocity  during  At.. 

J  J 


The  first  expression 


is  obtained  using  Equation  A-8.  It  is  arrived  at  by  replacing  x.,  by  r.  t. 

0  3  J 

and  solving  for  At..  The  second  expression  is  used  to  limit  the  fraction 

J 

of  the  foam  gasified  during  At ^  to  one-tenth. 

Burn  conditions  are  arrived  at  by  means  of  successive  approximations. 

Each  trial  r^  is  checked  using  the  resultant  values  for  the  foam  tem¬ 
perature  T^  and  the  heat  flux  q^  or  q  conducted  into  the  solid  propellant. 

The  trial  r.  value  is  checked  by  first  substituting  T^  into  Equation  8  to 
find  qo<  Then  the  boundary  condition  given  by  Equation  5  is  solved  for  r. 

If  r  does  not  agree  with  the  trial  value  r^  within  0.3  percent,  r^  is  revised 
and  the  calculations  repeated.  Usually  one  to  three  trials  are  needed  to 
achieve  the  above  accuracy. 

2.5.2  Gcu/Ouick  Pnoe.e.dun.t6 


the 


During  each  time  step  At.  the  gas/crack  predictions  involve  calculating 

J 

•  velocities  Wi,  W2  of  the  crack  walls 

•  crack  width  C 

w 

•  gas  density 

•  internal  energy  e 

•  gas  temperature  T^ 

•  gas  velocity 

•  gas  pressure  P 

•  * 

Velocities  Wi  and  Uz  of  the  crack  walls  are  computed  by  substituting  the 
bum  velocity  r^  and  the  pressure  P  into  Equations  26  and  27.  Changes  of  the 
crack  width  C  are  computed  by  substituting  the  velocities  U\  and  W2  into 

W 

Equation  28  and  multiplying  by  the  time  step. 
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3.  HMX  FOAM  MASSES 


3.1  FOAM  MASSES  DURING  STEADY  BURNING 

The  heat-transfer  constants  Ci  and  c2,  presented  in  Appendix  B  may  be 
used  in  conjunction  with  Equations  5  through  14  to  predict  the  foam  mass  during 
steady  burninq  as  a  function  of  the  pressure  P.  Resultant  foam  masses  are 
shown  in  Figure  4  in  terms  of  q/cm2.  Table  1  provides  tabular  results 
for  the  predicted  foam  masses  along  with  their  associated  steady  heat  fluxes 
qT-,  regression  rates  and  foam  temperatures  Tf.  Results  presented  in 
Figure  4  and  Table  1  are  predicted  on  neglible  lateral  gas  flows  across  the 
burning  propellant  surfaces. 

From  Figure  4  it  may  be  observed  that  the  foam  mass  decreases  with  in¬ 
creased  pressure.  The  latter  is  explained  by  the  consequence  of  elevated 
pressures  upon  the  rates  of  gasification  of  the  foam  and  upon  the  rate  of 
melting.  When  the  pressure  is  suddenly  Increased  to  a  fixed  higher  value, 
the  heating  of  the  foam  will  rise.  Increased  heating  causes  the  rate  of  gasi¬ 
fication  to  rise  more  rapidly  than  the  rate  of  melting  in  that  "excess  foam" 
accentuates  the  rate  of  heating  and  temperature  of  the  foam  directly  and  sub¬ 
sequently  the  melting.  Once  sufficient  form  is  expended,  the  rate  of  gasifi- 
catir,n  will  slow  while  the  entrance  of  relatively  cool  melt  acts  to  stabilize 
the  quantity  of  foam  at  a  lower  value. 

From  the  above  discussion,  it  should  be  clear  that  the  shape  as  well  as 
the  magnitude  of  the  curve  of  the  steady  foam  mass  versus  pressure  affects 
the  quantity  of  "excess  foam"  following  pressure  rises.  Rapid  evolution  of 
"excess  foam"  is  believed  to  be  the  single  most  important  factor  responsible 
for  the  production  of  pronounced  pressure  transients  needed  to  cause  detona¬ 
tion.  The  latter,  of  course,  remains  to  be  proven. 


I 


i 


I 


.  i 

y 


TABLE  1.  STEADY  HEAT  FLUXES  (qf),  BURN  RATES  (rf),  FORM  MASSES  (Mf), 
AND  FOAM  TEMPERATURES  (Tf)  DURING  HMX  BURNING  VERSUS  PRESSURE  (P) 


- - — . — - — -  —  — - - - 

P,  Bars 

qf,  cal/cm2-sec 

rf,  cm/sec 

g/cm2 

Tf,  K 

1 

10.0 

0.030 

0.0367 

601 

2 

19.0 

0.055 

0.0247 

615 

5 

44.0 

0.120 

0.0132 

635 

10 

84.0 

0.217 

0.0075 

653  j 

20 

160.0 

0.394 

0.0039 

673  j 

30 

233.0 

0.559 

0.0027 

686  1 

40 

305.0 

0.716 

0.0021 

695  1 

50 

376.0 

0.867 

0.0016 

702 

60 

446.0 

1.010 

0.0014 

709  j 

70 

515.0 

1.160 

0.0012 

714  j 

80 

583.0 

1.300 

0.0010 

719 

90 

652.0 

1.440 

0.0009 

723 

i 

100 

719.0 

1.530 

0.0008 

72/  1 

j 

1 

1 

> 

1 
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At  a  pressure  of  1  bar,  the  predicted  HMX  foam  mass  is  roughly  four 
times  greater  than  those  reported1  for  smaller  types  of  secondary  high  ex¬ 
plosives  such  as  Composition  B  and  PBX  9404.  Foam  masses  for  the  above  ex¬ 
plosive  materials  are  0.010  and  0.008  g/cm2,  respectively,  during  steady 
burning  at  1  bar  of  pressure. 

2.2  EXPERIMENTAL  VERIFICATION  OF  FOAM  MASS 

A  series  of  experiments  was  conducted  to  check  the  predicted  foam  mass 
of  burning  HMX  at  1  atmosphere  of  pressure.  These  experiments  used  a  HMX 
composite  propellant  described  in  Section  5.  The  foam  mass  was  determined 
by  measuring  the  quantity  of  sensible  heat  within  the  foam  layer,  and  divid¬ 
ing  the  result  by  the  product  of  the  burn  area  and  the  enthalpy  of  foam  during 
steady  burning.  The  foam's  heat  was  determined  b„  first  measuring  the  total 
sensible  heat  absorbed  by  burning  propellant  specimens,  and  then  substract- 
ing  the  experimentally  determined  heat  within  the  solid  portion  of  the  pro¬ 
pellant.  A  description  of  the  experiments  and  results  follows. 

3.2.1  H&U  Absorbed  by  BuAKtng 

Figure  5  illustrates  the  experimental  setup.  Burning  was  initiated 
by  igniting  a  thin  layer  of  ball  powder  sprindled  over  the  top  propellant 
surface.  An  insulated  steel  tube  was  used  to  contain  tho  resultant  flames. 
Thermocouples  (5  mil  chromel-alumel )  were  imbedded  at  various  depths  in  the 
upper  half  of  the  propellant.  These  thermocouples  served  to  measure  the  rate 
of  burning,  and  determine  the  temperature  profile  within  the  solid  propellant. 
After  steady  burning  was  achieved,  the  remainder  of  the  burning  propellant 
disk  was  dropped  into  the  water  bath  contained  within  the  calorimeter.  Burn¬ 
ing  ceased  immediately  upon  entry  of  the  disk  into  the  water.  Then  the  calo¬ 
rimeter  was  closed  with  a  cork  stopper  housing  a  Berkman  Thermometer.  This 
thermometer  measured  the  resultant  temperature  rise  to  within  ±0.01°C. 

Next  the  total  head  capacity  Hc  of  the  water  bath,  calorimeter  and  ther¬ 
mometer  was  measured  by  adding  a  known  quantity  of  heat  to  the  calorimeter, 
and  observing  the  resultant  temperature  rise.  In  each  test,  the  initial  tem¬ 
peratures  of  all  media  wore  within  0.5°C  of  that  of  the  environment. 
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The  total  heat  Q  absorbed  per  unit  area  of  burn  propellant  was  determined 


Q  »  [(He  +  M  Cp)  AT  -  Qf]/ A  (30) 

where: 

Hc  *  measured  heat  capacity  of  calorimeter  plus  water,  cal/°C 

M  *  mass  of  propellant,  g 

Cp  =  specific  heat  of  propellant,  cal/gm-°C 

AT  =  temperature  tise  of  calorimeter,  °C 

Qf  =  heat  transferred  to  calorimeter  during  free-fall  of 
burning  propellant,  cal 

A  *  area  of  propellant’s  burn  surface,  cm2 

The  specific  heat  Cp  of  the  propellant  was  measured  at  a  temperature  of 
24°C.  It  was  0.25  cal/g-°K.  The  quantity  of  heat  Qf  transferred  from  flames 
to  the  calorimeter  was  found  by  dropping  another  burning  propellant  disk  into 
an  insulated  vessel  containing  water  which  was  partially  submerged  within  the 
water  bath  of  the  calorimeter.  The  latter  vessel  was  withdrawn  along  with  the 
propellant's  sensible  test  immediately  after  burninq  was  extinguished  by  the 
water.  Then  the  calorimeter  was  closed  before  observing  the  temperature  rise 
of  the  calorimeter  caused  by  the  propellant  flames.  This  heat  equalled  6  per¬ 
cent  of  the  total  heat  absorbed  by  the  calorimeter  when  the  propellant's  heat 
was  not  withdrawn  from  the  calorimeter. 

The  total  heat  Q  absorbed  by  the  propellant  disk  during  steady  burnlno 
at  1  atmosphere  of  pressure  was  12.4  cal  per  cm2  of  burn  surface.  The  latter 
represents  the  mean  of  3  measurements  given  by  10.5,  11.7  and  15.1  cal/cm2. 
These  measurements  include  sensible  heats  within  the  foam  and  underlying 
solid  propellant.  The  burning  rate  was  0.3  cm/sec. 


3.2.2  Heai  in  Solid.  PKopzJULant  Vu 


8 u/inir 


The  sensible  heat  in  the  solid  propellant,  which  is  represented  by  <5Q, 
was  determined  from  the  thermocouple  measurements  mentioned  earlier.  The 
thermocouples  provided  continuous  temperature-time  curves  within  the  solid 
portion  of  the  burning  propellant  as  the  burning  approached  steady  velocities. 
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Only  those  portions  of  the  temperature-time  curves  associated  with  steady 
burning  were  used.  The  steady  temperature  profile  T(x)  versus  propellant 
depth  x  was  determined  using  the  above  temperature-time  curve  in  conjunction 
with  the  rate  of  burning.  The  heat  was  determined  In  terms  of  the  enthalpies 
H(T)  corresponding  to  these  temperatures  T(x)  as  follows: 

oo 

<$Q  »  pj~  H(T(x))  dx  (31) 

o 

The  values  of  H(T)  used  in  equation  31  are  those  for  HMX  shown  In  Figure  6. 
They  were  obtained  from  Reference  1.  -5Q  was  evaluated  as  4.3  cal  per  cm 

of  burn  surface. 

3.2.3  Steady  Foam  Masa  cut  1  ktmoApheJie.  o &  PKtUiuAe. 

The  foam  mass  was  determined  as  follows: 


The  foam  temperature  Tf  cited  above  equals  that  to  support  the  steady  burning 
rate  of  HMX  at  1  atmosphere.  The  latter  temperature  was  predicated  as  6lO°K 
from  Table  1.  The  Initial  temperature  Tn  was  294  K.  The  enthalpies  H  are 
assumed  equal  to  those  of  HMX  shown  in  Figure  6. 

The  three  mass  determinations  were  0.035,  0.042,  and  0.062  a/cm2.  The 
mean  of  the  three  values,  namely  0.046  g/cm2  is  in  rough  agreement  with  the 
value  0.037  presented  by  Figure  4  for  pure  HMX  at  1  bar  of  pressure.  Thus 
the  curve  of  Figure  4  roughly  approximates  three  experimental  points  at  1, 

34  and  68  bars. 

Whether  foam  exists  at  pressures  above  68  bars  or  sublimes  away  remains 
In  question.  In  this  regard  foam  can  be  generated  within  burning  propellant 
cracks  ut  pressures  between  1  and  68  bars.  Analysis  presented  in  Section  4 
Indicates  that  rapid  gasification  of  such  foam  can  cause  pronounced  pressure 
rises  within  cracks;  and  lead  to  DDT  provided  sufficient  foam  Is  generated. 
Rapid  sublimation  of  the  foam  at  pressures  above  68  bars  would  accentuate 
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Figure  6.  Enthalpy  of  HMX  versus  temperature. 


the  pressure  rises  and  make  DDT  more  likely.  Therefore,  In  assuming  the 
existence  of  foam  at  all  pressures,  we  believe  we  are  conservative  insofar 
as  DOT  Is  concerned. 
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4.  MODEL  PREDICTIONS 


Before  discussing  model  predictions  it  Is  important  to  appreciate  the 
magnitudes,  and  durations  of  short  wave  pressures  needed  to  cause  DDT.  In 
this  regard  de  Longuiville6  presents  threshold  shock  wave  amplitudes  (P) 
and  durations  (t)  associated  with  the  initiation  of  HMX  -nylon  by  Impact.  His 
results  are  presented  in  Figure  7.  The  curve,  which  represents  a  constant 
weighted  impulse  P2t,  separates  conditions  in  which  detonation  did  and  did 
not  occur.  Clearly  weighted  impulse  has  limited  usefulness  in  that  it  ap¬ 
plies  to  step-wise  pressure  waves  and  durations  cited  in  the  figure.  Never¬ 
theless  de  Longuiville' s  results  suggest  DDT  will  require  pressure  rises  of 
the  order  of  tens  of  K  bars  or  more  within  times  of  the  order  of  a  y  sec  or 
less.  The  more  abrupt  the  pressure  rise  is  and  the  slower  its  decay,  the 
greater  the  likelihood  of  DDT  is. 

4.1  SENSITIVITY  ANALYSES 

Here  we  are  concerned  with  examininq  the  effect  of  various  factors  upon 
pressure  transients  produced  in  burning  cracks.  First  the  consequence  of  al¬ 
tering  various  propellant  properties  and  crack  conditions  upon  pressure  tran¬ 
sients  within  single  cracks  is  examined.  Then  the  consequence  of  applying 
pressure  transients  (or  stress  waves)  from  one  crack  to  the  next  in  a  sequen¬ 
tial  fashion  is  determined.  Table  F-l  in  the  Appendix  F  describes  property 
values  used  for  the  HMX  propellant  and  combustion  gases. 

Six  parameters  were  varied  in  this  study.  Two  of  the  six  parameters  are 
propellant  properties.  These  are: 

•  internal  heat  Qs 

•  propellant  impedance  I0  at  ambient  pressure 
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Figure  7 .  Impact  results  from  Reference  6  for  HMX-nylon  explosion 


] 

4 

a  o 

i 

»  . 

j. 

i 

t 


L 


Remaining  parameters  are: 

•  initial  crack  width  Cwo 

•  initial  gas  pressure  P0  (foam  temperature 
adjusted  accordingly) 

•  amplitude  AP  of  incident  stress  wave 

•  initial  foam  mass  Mf0 

The  adjective  "initial"  refers  to  values  immediately  before  the  stress  wave 
arrives.  Each  of  the  latter  parameters  can  vary  widely  from  crack  to  crack 
depending  upon  how  the  crack  develops,  ignites  and  burns.  Uncertainties  also 
exist  in  the  propellant  properties.  For  these  reasons,  three  values  were 
chosen  for  each  of  the  six  parameters  cited  above.  They  are  listed  in  Table 
2.  In  each  study  the  initial  foam  temperature  was  set  equal  to  its  value 
during  steady  burning  at  the  initial  pressure  P0. 
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Values  of  the  internal  heat  Q$  and  propellant  Impedance  IQ  were  estimated. 
The  Qs  values  are  slightly  smaller  than  cited  in  the  literature7  for  HMX. 

The  values  for  the  initial  gas  pressure  PQ  and  the  amplitude  AP  of  the  inci¬ 
dent  stress  wave  are  considered  typical  of  the  ma3iutudes  one  may  encounter 
in  a  crack.  In  this  regard  PQ  varies  with  time  and  gas  flow  into  a  crack 
while  AP  depends  upon  the  cavity  pressure  and  rocket  motor  geometries. 

Foam  masses  presented  in  Table  2  are  larger  than  the  values  present¬ 
ed  In  Figure  4  for  steady  burning  at  the  pressure  P  .  Implicit  in  this  as¬ 
sumption  is  relatively  low  heating  rates  prior  to  ignition,  or  the  cumulation 
of  melt  due  to  melt  flow  brought  about  by  gas  flows  in  the  crack. 

Figure  8  is  presented  to  illustrate  salient  features  of  the  problem.  It 
presents  transient  crack  widths,  melt  passes  and  pressures  using  the  unper¬ 
turbed  values  presented  in  Table  2.  Time  starts  with  the  arrival  of  the  given 
stress  wave. 

Initially  the  stress  wave  partially  collapses  the  crack.  The  result  ir 
increased  gas  pressures  which  cause  increased  propellant  heating  (see  Equations 
11  and  13).  The  result  is  accelerated  burning  that  supports  progressive  in¬ 
creases  of  the  pressure. 

Early  during  the  pressure  buildup  the  crack  commences  to  expand  in  res¬ 
ponse  to  the  elevated  pressures.  Pressures  continue  to  rise  during  crack  ex¬ 
pansion  until  the  "excess  melt"  is  consumed  by  burning.  The  pressure  spikes 
are  due  to  the  fact  that  much  of  the  melt  is  consumed  within  short  times  of 
the  order  of  a  few  p  seconds.  Thereafter  the  burning  rates  and  pressures 
comnence  to  decrease  due  to  continued  crack  expansion. 

Figure  9  indicates  that  higher  internal  heats  Qs  promote  higher  pressures 
at  earlier  times.  The  latter  is  due  to  greater  propellant  heating  with  higher 
Qg  values.  Figure  10  shows  that  higher  propellant  Impedances  IQ  also  cause 
higher  pressures.  In  this  case  higher  pressures  are  due  to  decreases  In  crack 
expansion.  .  t 

In  Figure  11  the  initial  pressure  PQ  and  foam  temperature  were  varied 
simultaneously.  It  shows  that  increased  initial  pressures  and  foam 
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Figure  9.  Consequences  of  varying  internal  heat  ( Q - 
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Figure  11. 
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temperatures  result  in  greater  pressures  at  earlier  times.  Increased  pres¬ 
sures  are  due  to  reductions  in  crack  expansion  caused  by  more  rapid  consump¬ 
tion  of  the  "excess  foam". 


Smaller  crack  volumes  can,  of  course,  also  be  achieved  by  the  amplica¬ 
tion  of  higher  amplitude  stress  waves  or  by  starting  with  smaller  cracks. 
Their  effects  are  shown  in  Figures  12  and  13.  Of  the  two,  stress  waves  are 
considered  more  important  in  that  it  is  possible  to  generate  stress  waves 
orders  of  magnitude  greater  than  those  cited  in  Figure  12.  In  addition  it  is 
difficult  to  maintain  extremely  small  crack  widths  unless  such  stress  waves 
are  applied  to  counter  crack  expansion  by  the  pressure  buildup.  Figure  13 
shows  that  smalle-  cracks  grow  at  significantly  greater  rate  than  larger 
cracks. 

Figure  14  shows  that  larger  initial  foam  masses  produce  pronounced 
increases  of  the  pressures.  There  are  two  reasons  for  the  higher  pressures. 
The  first  is  the  larger  amounts  of  gas  evolved.  The  second  is  the  more  rapid 
gasification  of  the  foam  layer.  The  latter  is  explained  by  the  cooling  ef¬ 
fects  of  the  molten  propellant  entering  the  foam.  Temperature  rises  are  in¬ 
hibited  less  by  the  incoming  melt  with  larger  foam  masses.  The  result  is 
more  rapid  gasification  of  the  "excess  foam".  In  turn  rapid  gasification 
yields  less  time  for  crack  expansion  and  hence  higher  pressures. 

Next  let  us  select  three  sets  of  values  for  the  six  parameters  presented 
in  Table  2  to  gain  a  better  appreciation  of  the  range  of  pressures  that  may 
be  produced  ir,  cracks.  The  three  sets  of  values  are  presented  in  Table  3. 

The  case  2  values  represent  nominal  values  presented  in  Table  2.  The 
case  1  values  are  those  that  yielded  the  smallest  pressures;  the  case  3 
values  are  those  that  yielded  the  highest  pressures. 
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TABLE  3.  PARAMETRIC  VALUES  FOR  THREE  SELECTED  CRACK  CONDITIONS 


Parameter 

Case  1 

Values 

Case  2 

Values 

Case  3 

Values 

Qs 

100  cal/g 

150  cal/g 

200  cal/g 

*0* 

0.30  bars  sec/cm 

0.45  bars  sec/cm 

0.60  bars  sec /cm 

^wo 

0.15  cm 

0.10  cm 

0.05  cm 

AP 

34  bars 

68  bars 

102  bars 

Po+ 

17  bars 

34  bars 

51  bars 

Mfo 

0.005  q/cm2 

0.010  g/cm2 

0.015  g/cm2 

*  Impedance  I  =  IQ  (1+0.0002  P(bars)) 

f  Foam  temperature  and  temperature  distribution  of  solid  propellant 
varied  with  PQ  accordino  to  steady  burning  conditions 

Resultant  pressure  transients  are  presented  in  Figure  15  for  the  three 
cases.  Notice  that  peak  pressures  differ  by  an  order  of  magnitude.  It  sug¬ 
gests  that  pressure  transients  will  vary  widely  from  crack  to  crack  depending 
upon  the  propellant  properties  and  how  the  crack  develops,  ignites  and  burns. 

4.2  MULTIPLE  CRACKS 

Multiple  burning  cracks  are  important  in  that  the  pressures  produced 
by  stress  wave/crack  interactions  are  greater  than  the  applied  stress  wave. 

The  result  is  an  enhanced  stress  wave  leaving  cracks  that  can  then  act  upon 
neighboring  burning  cracks.  By  this  process  prcgressly  higher  and  steeper 
pressure  transients  can  develop  as  stress  waves  more  from  crack  to  crack. 

To  better  appreciate  the  problem  consider  three  sets  of  identical  parallel 
cracks.  The  cracks  are  identical  to  those  considered  earlier  in  the  sec¬ 
tion.  They  are  described  by  case  1,  case  2,  and  case  3  in  Table  3. 
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Time  Following  Arrival  of  Stress  Wave 
.fW*  IS.  Pressure  transients  for  three  crack  conditions 


Figure  16  Illustrates  the  results  for  a  series  of  Identical  case  1  cracks. 
Notice  that  the  pressure  transients  become  more  pronounced  with  each  succed- 
Ing  crack.  Much  of  the  pressure  Increase  is  caused  by  progressive  decreases 
of  the  crack  width.  At  least  seven  cracks  are  needed  to  achieve  pressures 
of  the  order  of  10  k  bars. 

Figures  17  and  18  present  similar  results  for  case  2  and  3  cracks,  res¬ 
pectively.  In  each  of  the  above  cracks,  pressures  rise  more  steeply  to  their 
peaks  than  the  pressures  produced  by  the  case  1  cracks  (see  Figure  16). 

Steep  pressure  rises,  are  of  course,  more  conductive  to  detonation.  Also 
fewer  cracks  are  needed  to  develop  pressures  of  the  order  of  tens  of  kbars 
with  case  2  and  3  cracks  than  required  with  case  1  cracks. 

The  results  of  Figures  16,  17,  and  18  suggest  that  multiple  cracks  may 
lead  to  detonation  provided  enough  cracks  are  involved  in  the  sequential 
fashion  indicated.  This  hypothesis  suggest  that  large  propellant  motors  are 
more  susceptible  to  detonation  than  small  motors.  It  is  consistent  with  un¬ 
reported  IITRI  observations  in  which  a  few  hundred  pounds  of  secondary  high 
explosives  (HE)  burned  freely  without  event,  while  a  few  thousand  pounds  of 
the  same  HE  detonated  under  similar  burning  conditions.  Each  test  result 
was  replicated  four  times. 

4.3  FOAM  MASSES  GENERATED  WITHIN  CRACKS 

Propellant  surfaces  in  cracks  are  exposed  to  cross  gas  flows  that  accen¬ 
tuate  propellant  heating  and  therby  lessen  foam  formation.  In  this  section 
we  shall  present  prodicted  foam  masses  generated  in  cracks  that  propogate 
into  a  cavity  of  high-temperature  high-pressure  gases.  Initially  the  gases 
in  the  cracks  are  assumed  to  have  a  temperature  of  294°K  and  a  pressure  of 
1  bar.  The  downstream  end  of  the  crack  is  considered  connected  to  a  second 
cavity  of  gases  held  at  the  same  temperature  and  pressure  as  the  gases  ini¬ 
tially  in  the  crack. 
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Figure  16.  Multiple  cracks  lease  1)  concluded. 
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Figures  19  through  21  present  dynamic  foam  masses  as  well  as  gas  velo- 
cities,  cracks  widths  and  gas  pressures  along  the  length  of  the  crack  for 
cavity  pressures  of  8,  16  and  68  bars,  respectively.  Conditions  are  pre¬ 
sented  for"three  times  following  exposure  of  the  crack  to  the  high-pressure 
cavity. 

As  one  would  expect,  the  foam  mass  is  greatest  initially  at  the  high- 
pressure  end  of  the  crack.  Thereafter,  the  qreatest  foam  mass  occurs  at 
ncreasing  crack  depths  before  being  found  at  the  downstream  end  of  the 
crack.  The  latter  is  due  to  the  smaller  downstream  heat  fluxes; 
i.e.,  the  heat  fluxes  at  the  upstream  end  are  relatively  high.  In  this 
regard  lower  fluxes  produce  greater  foam  masses  over  longer  periods  of  time. 

Mean  values  of  the  foam  masses  within  the  crack  are  presented  in  figure 
22  as  a  function  of  time  for  each  of  the  three  cases  illustrated  in  Figures 
19  through  21.  It  may  be  observed  that  the  mean  foam  mass  asymototically  ap¬ 
proaches  values  of  approximately  0.0040,  0.0020,  and  0.0005  q/cm2  for  cavity 
pressures  of  8,  16,  and  68  bars,  respectively.  These  foam  masses  are  roughly 
40  percent  of  their  steady-state  values  presented  in  Figure  4.  Figure  22 
indicates  that  greater  amounts  of  foam  require  longer  times  to  produce. 

The  relatively  low  foam  masses  cited  above  do  not  appear  adequate  to 
cause  DDT  unless  extremely  intense  stress  waves  are  generated  say  by  multi¬ 
ple  burning  cracks.  Two  possibilities  exist  for  foam  enhancement.  The 
first  and  most  obvious  is  low  gas  pressures  and  propellant  heating  rates, 
this  possibility  does  not  appear  likely.  The  second  possibility  is  accumu¬ 
lation  of  melt  at  particular  regions  of  a  crack  due  to  flow  of  melt  caused 
by  high  velocity  gases.  This  possibility  is  more  likely  in  that  mixing 
of  relatively  cool  melt  with  hotter  melt  will  lessen  the  net  gasification  of 
the  melt,  and  thereby  increase  the  total  amount  of  melt  present. 
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Figure  21.  Consequence  of  sudden  exposure  of  "open  crack"  to  high 
pressure  (68  bars )  cavity. 
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5.  BURNING  ,'ROPELLANT  EXPERIMENTS 


5.1  EXPERIMENTAL  PROCEDURE 

A  series  of  18  experiments  were  conducted  with  the  objective  of  validating 
the  analytical  conclusions  described  in  Section  4.  A  schematic  illustrating 
essential  features  of  the  experimental  set-up  is  presented  in  Figure  23.  The 
design  simulates  partial  closure  of  burning- propellant  cracks  by  stress  waves. 
Nomenclature  is  included  in  Figure  23  for  future  reference. 

In  the  experimental  set-up  of  Figure  23,  a  lead  driver  is  dropped  onto 

the  piston  assembly  which  closes  the  vents  as  it  moves  downward.  Pressure 
then  rises  due  to  reduction  of  the  void  space  and  increases  of  the  burning 
rate  brought  about  by  increased  pressure.  A  necessary  but  not  sufficient 
condition  for  DDT  is  reductions  of  the  height  of  the  void  space  to  hundreths 
of  a  cm  or  less  in  order  to  generate  pressures  of  the  order  of  ten  k  bars  or 
more. 

Ignition  was  accomplished  by  placing  six  0.5  cm  diameter  balls  of  the 
propellant  along  with  4  gm  of  ball  powder  upon  the  upper  surface  of  the  pro¬ 
pellant  cylinder  located  at  the  bottom  of  the  chamber.  The  ball  powder  was 

ignited  by  an  electrical  heater  wire.  Then  the  burning  ball  powder  ignited 

the  balls  of  propellant  that  ignited  the  cylinder  of  propellant.  In  this 
regard,  ball  powder  bums  too  rapidly  to  ignite  the  propellant  cylinder 
directly.  Propellant  balls  are  more  susceptible  to  ignition  than  the 
surface  of  the  propellant  cylinder  because  much  of  their  surface  is  in  the 
flame  zone.  Immediately  following  ignition,  the  electrical  heater  wire  was 
withdrawn  from  the  chamber  so  it  would  not  interfere  with  subsequent  move¬ 
ments  of  the  piston. 

Following  igntion,  the  propellant  was  allowed  to  burn  for  25  sec  to 
achieve  steady  burning.  Steady  burning  is  desired  in  that  it  is  easier  to 
characterize  than  transient  burning.  Then  the  driver  was  released  forcing 
the  piston  assembly  into  the  chamber.  The  nominal  clearance  between  the 
bearing  surfaces  of  the  piston  and  chamber  was  only  0.001  in.  to  deter  gas 
escape.  Various  lubricants  were  used  to  maintain  piston  velocities. 
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Figure  27.  Schematic  of  experimental  (setup. 
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Closest  approach  of  the  burning  propellant  surfaces  (or  minimum  height 
of  the  void  space)  was  determined  by  attaching  a  metal  pin  to  the  lower  face 
of  the  piston.  This  pin  extended  through  any  propellant  disk  attached  to  the 
piston  so  the  pin  would  indent  the  propellant  cylinder.  Such  indentations 
would  be  preserved  by  sudden  extinguishment  of  the  burning  provided  the  explo¬ 
sion  did  not  destroy  the  propellant  cylinder.  Extinguishment,  of  course,  is 
caused  by  reduction  in  the  propellant  heating  brought  about  by  sudden  pressure 
relief,  and  by  inadequate  heat  stored  in  the  propellant  to  maintain  burning 
at  ambient  pressure. 

In  the  remainder  of  this  section,  we  shall  describe  the  experiments  and 
analysis  used  to  interprete  the  experimental  results. 


5.2  ANALYSIS 

Three  important  factors  affect  the  response  of  the  burning  propellant. 
The  first  is  increases  of  the  heat  fluxes  brought  about  by  pressure  rises. 
This  aspect  of  the  problem  is  described  in  Section  2.1.1.  The  second  factor 
is  velocities  of  the  piston  assembly  and  driver  following  impact.  Impacts 
were  assumed  ideal  in  which  energy  and  mementum  is  conserved.  Resultant 
velocities  Vp  and  of  the  piston  assembly  and  ariver,  respectively,  are 


V 

V 


P 

d 


(2  Md  vdo  4  (Mp  •  V  V/!Md  4  V 


»  HP  vpo  4  <"d  -  V  V'^d 4  v 


(33) 

(34) 


where  Mp  and  represent  masses  of  the  piston  assembly  and  driver,  respec¬ 
tively,  while  Vp0  and  VdQ  represent  thier  velocities  immediately  prior  to 
impact.  The  above  equations  not  only  apply  to  the  initial  impact  wherein 
Vpo  =  0,  but  also  to  any  subsequent  impacts. 

The  third  and  last  feature  of  the  experiments  has  to  do  with  drag  forces 
imparted  the  moving  piston  by  lubricants  and  contaminant  produced  by  the 
burning.  These  drag  forces  are  assumed  to  vary  linearly  with  the  piston 
velocity  so  that 
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where  and  Ap  represent  the  areas  of  the  bearing  surface  and  face  of  the  piston, 
respectively,  while  d  and  g  represent  the  drag  coefficient  and  acceleration 
of  gravity,  respectively.  Pressure  rises  per  unit  area  are  represented  by 
AP.  Preliminary  estimates  of  the  drag  coefficients  d  were  determined  prior 
to  the  burning  propellant  experiments  by  measuring  the  terminal  velocity  Vp 
of  the  piston  and  applying  the  applicable  terms  of  Equation  35  to  yield 

d=Mpg/(AdVp)  (36) 


Actual  d  values  were  estimates  by  use  of  a  computer  code  that  accounts  for 
each  of  the  above  phenomena.  Evaluation  was  made  by  trial  and  errur  until 
predicted  closure  of  the  void  space  agreed  with  measurements.  Drag  co¬ 
efficients  were  along  appreciably  higher  than  preliminary  estimates  due  to 
contamination  and/or  loss  of  lubricant. 


5.3  EXPERIMENTAL  RESULTS 

A  composite  propellant  containing  appreciable  HMX  was  selected  due  to 
the  lack  of  apparatus  to  press  pure  HMX  cylinders  of  the  large  sizes  desired. 
The  propellant  samples  were  manufactured  by  the  Jet  Propulsion  Laboratory 
(AFOPL)  in  California.  The  AFJPL  specifications  indicated  that  the  composite 
propellant  contained: 

PEG  6.250  percent 

CAB  0.250 

N-100  12.500 

TEG0N  12.500 

HMXa  40.000 

HMXe  25.000 

MNA  1.000 

ZRC  1.000 

CARBON  0.500 

TPB  0.025 

MALIC, 

ANHY.  0.025 _ 

*100.600  percent 

According  to  AFJPL,  the  resultant  propellant  cylinders  were  of  poor  quality. 

They  were  very  oily,  stuck  to  anything  they  touched,  and  did  not  hold  their 
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shape  when  removed  from  their  molds.  In  order  to  achieve  firm  propellant 
cylinders,  the  propellant  was  wrapped  in  highly  absorbent  paper  to  remove 
the  excess  oil.  Then  the  propellant  was  pressed  into  desired  shapes  to 
eliminate  any  velds.  In  doing  so,  the  propellant  specimens  incurred  a  weight 
loss  of  approximately  7  percent. 

Eighteen  experiments  were  conducted  using  the  experimental  set-up 
illustrated  by  Figure  23.  The  first  two  of  the  eighteen  experiments,  namely 
Nos.  1  and  2,  were  used  to  develop  the  ignition  system  described  earlier  in 
Section  5.1.  Table  4  summarizes  the  remaining  experiments.  Nomenclature  is 
presented  in  Figure  23.  The  experiments  differed  according  to  the: 

•  height  and  weight  of  the  driver 

•  weight  of  piston  assembly 

•  presence  of  propellant  disk  on  piston 

•  location  of  propellant  surfaces  with  respect  to 
vent  holes 

•  type  of  lubricant  used  on  bearing  surfaces 

In  tests  3  through  8  the  bearing  surfaces  were  coated  with  a  0.002  in. 
film  of  silicone  grease  (prior  to  inserting  piston  into  chamber),  and  the  driver 
weights  and  heights  were  progressively  increased  in  an  attempt  to  increase  the 
order  of  the  explosion  by  improving  closure  of  the  void  space.  Test  4  was 
lost  due  to  failure  of  the  driver  release  mechanism.  Tests  6  and  7  produced 
the  most  severe  explosion  driving  the  101  lb  driver/piston  against  the  rein¬ 
forced  concrete  ceiling  of  the  test  facility.  The  explosion  was  clearly  not 
a  detonation.  In  fact  the  only  damage  to  the  chamber  and  piston  was  caused 
by  the  falling  driver  and  piston  assembly.  Test  8  yielded  the  least  severe 
explosion  of  the  6  tests.  In  fact  there  were  two  minor  explosions.  The  first 
explosion  was  due  to  venting  caused  by  the  piston  being  driven  above  the  vents 
by  the  pressure  rise.  The  second  explosion  was  caused  by  the  driver  driving 
the  piston  downward  pass  the  vent  holes  a  second  time. 

The  greatest  closure  or  smallest  height  of  the  void  space  obtained  in 
tests  3  through  8  was  0.8  cm  in  test  7.  All  other  closures  were  in  excess  of 
1.0  cm.  Such  void  space  heights  are  much  too  great  to  achieve  DDT.  Initially, 
the  lack  of  closure  was  attributed  solely  to  burning  propellant  particles 
becoming  wedged  between  the  piston  and  chamber  after  being  embedded  in  the 
lubricant.  ||T  #£SEAICH  institute 
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As  a  consequence,  it  was  decided  to  conduct  test  9  with  no  lubricant 
whatsoever.  It  may  be  observed  that  closure  was  not  improved.  In  fact,  Sx 
Increased  to  1.5  cm.  Residues  found  on  the  chamber  walls  after  the  test  pre¬ 
vented  a  clean  piston  from  being  pushed  through  the  chamber  by  hand. 

At  this  point,  it  was  decided  to  sprinkle  powdered  graphite  upon  the 
silicone  grease  film  used  in  tests  3  through  8.  It  was  hoped  that  the  powdered 
graphite  would  deter  combustion  products  from  sticking  to  the  silicone  grease. 
This  lubricant  was  used  in  tests  10  and  11  without  success.  Drag  coefficients 
d  were  essentially  the  same  as  those  obtained  using  the  silicone  grease  without 
powdered  graphite.  Improved  closure  (0.6  and  0.7  cm)  of  the  void  space  was 
attributed  to  increased  driver  weight  and  height  along  with  a  greater  piston 
assembly  mass.  Even  though  the  closure  was  about  an  order  of  magnitude  less 
than  desired,  the  explosions  of  test  10  was  sufficiently  intense  to  propel 
209  lbs  cgainst  the  facility's  concrete  ceiling  with  considerable  force.  The 
guide  rod  shown  in  Figure  23  was  driven  into  the  concrete  to  a  depth  of  1  inch. 
Remainder  of  the  rod  was  bent  much  like  a  pretzel.  Energy  needed  to  propel 
the  piston/dri ver  against  the  ceiling  exceeded  35  percent  of  the  energy  stored 
In  the  foam  layer. 

In  tests  12  and  13,  it  was  decided  to  discard  the  use  of  silicone  grease 
in  favor  of  powdered  graphite  wetted  with  silicone  oil  for  purposes  of  In¬ 
tegrity.  Drag  coefficients  were  higher  than  those  found  earlier;  explosions 
were  less  severe  and  closure  was  not  significantly  Improved  from  tests  10  and 
11.  The  latter  is  inspite  of  greater  driver  heights.  It  is  suspected  that 
much  of  the  lubricant  was  swept  away  by  the  flow  of  combustion  products. 

The  best  lubricant  found  involved  using  an  extremely  thin  layer  of 
silicone  grease  to  hold  the  graphite  particles.  The  thin  silicone  grease 
film  was  formed  by  first  applying  the  grease  to  the  piston  and  chamber  walls, 
and  then  inserting  the  piston  into  the  chamber  to  remove  excess  grease.  This 
process  was  repeated  until  the  piston  fell  freely  through  the  chamber.  Then 
powdered  graphite  was  sprinkled  upon  those  lubricant  surfaces  that  subsequently 
will  be  exposed  to  the  hot  combustion  products.  This  lubricant  was  first  used 
in  test  14.  Two  reasons  may  be  advanced  for  the  improved  lubrication.  These 
are  greater  resistance  of  the  film  to  the  flow  of  combustion  gases,  and  a 
reduction  of  the  liklihood  of  particles  being  embedded  In  the  thin  film  of 
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silicone  grease.  These  reaons  are  of  course,  conjectural.  The  above  lubricant 
was  used  In  each  of  the  remaining  tests. 

In  test  15,  the  masses  of  the  driver  and  piston  assembly  were  increased 
appreciably,  and  the  height  Sx  of  the  void  space  (distance  between  bottom  of 
vent  holes  to  surface  of  burning  propellant  cylinder  at  time  of  driver  release) 
was  reduced.  This  test  was  particularly  discouraging  in  that  there  was  no 
Improvement  in  closure  and  the  driver/piston  assembly  was  thrown  less  than 
1  ft  upward. 

Two  possibilities  remained  for  the  inadequate  closure.  The  first  is 
thermal  distortion  of  the  exposed  chamber  walls;  and  the  second  is  much  more 
rapid  consumption  of  the  foam  layer  than  anticipated  by  analysis.  Test  16 
was  used  to  check  the  former  hypothesis.  The  sole  difference  with  test  15  was 
a  0.020  in.  reduction  of  the  piston's  diameter  so  the  piston  would  hopefully 
not  be  slowed  or  stopped  by  distortions  of  the  heated  chamber  walls.  In  this 
regard,  the  exposed  chamber  walls  absorbed  an  average  heat  flux  of  4  cal/cm2- 
sec,  based  upon  thermal  measurements.  The  essentially  identical  results  of 
tests  15  and  16  dicounts  thermal  distortions  of  the  chamber  walls  as  the 
prime  cause  of  poor  closure.  A  further  test  of  the  hypothesis  was  achieved 
by  redesigning  the  set-up  of  Figure  23  as  illustrated  by  Figure  24.  In  this 
set-up,  the  chamber  walls  are  thermally  protected  by  attaching  a  steel  sleeve 
to  the  piston  with  the  propellant  contained  within  the  sleeve.  The  clearance 
between  the  sleeve  and  piston  was  Increased  from  its  normal  value  of  0.001  in. 
to  0.004  in.  to  counter  any  thermal  expansion  of  the  sleeve.  This  device  was 
used  In  test  17.  The  lack  of  Improved  closure  suggests  thermally- induced 
distortion  was  not  the  cause  of  piston  stoppage. 

At  this  point,  it  concluded  that  lack  of  closure  was  due  to  the  production 
of  higher  than  anticipated  gas  pressures  while  the  piston  was  moving  downward. 
Two  possible  causes  are  ejection  of  molten  HMX  and  particulates  from  the  foam 
layer  into  the  combustion  zone,  or  a  more  rapid  increase  in  the  rate  of  Inter¬ 
nal  heating  than  predicted  by  pQ$rf.  Both  would  speed  gasification  of  the 
foam  layer  and  hence  accelerate  pressure  rises.  Unfortunately,  elaborate 
Instrumentation  needed  to  test  this  hypothesis  were  not  Included  in  the  scope 
of  work  and  could  not  afforded  by  existing  funds.  Instead,  it  was  decided 
to  vent  some  of  the  gas  during  closure,  and  determine--whether  or  not  closure 
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Is  improved  sufficient  to  cause  DDT.  The  additional  vent  was  provided  in  the 
final  test,  namely  No.  18.  An  illustration  indicating  provisions  for  the 
additional  venting  is  presented  in  Figure  25.  Here  a  single  0.25  in.  diameter 
vent  hole  was  located  beneath  the  surface  of  the  propellant  cylind..  with  a 
passage  way  cut  in  the  propellant.  Vent  size  was  purely  a  guess  lacking  in¬ 
formation  regarding  how  much  gas  must  be  allowed  to  escape  to  achieve  desired 
closures  of  the  order  of  hundreths  cf  a  cm  or  less. 

Remaining  features  of  test  18  are  summarized  in  Table  4  along  with  experi 
mental  results.  This  test  produced  two  low-order  explosions  very  similar  to 
that  produced  in  test  8.  The  composite  weight  of  275  lbs  of  the  driver  and 
piston  assembly  rose  less  than  1  ft.  Closure,  however,  improved  from  a  pre¬ 
vious  low  of  0.6  cm  to  0.3  cm.  Unfortunately,  the  resultant  closure  was 
roughly  an  order  of  magnitude  larger  than  predicted  fcr  DDT  by  analysis. 
Nevertheless,  the  improved  closure  was  sufficient  to  indicate  that  the  foam 
mass  is  consumed  much  faster  than  expected.  At  this  point  in  time,  further 
testing  was  suspended  due  to  the  lack  of  propellant. 
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6.  SUMMARY / CONCLUS IONS 


6.1  ANALYSES 

This  study  concludes  that  burning  HMX  melts  ana  forms  a  melt  on  a  foam 
layer  over  the  burning  surface.  During  steady  burning,  this  layer  progres¬ 
sively  '’'"'creases  with  increased  pressure.  The  primary  threat  of  DDT  of  HMX 
based  '■ellants  lies  in  the  relatively  high  thermal  energy  content  of  its 
melt.  When  triggered  into  dynamic  burning,  the  melt  will  rapidly  gasify.  In 
restricted  spaces  such  as  cracks,  its  gasification  will  generate  pronounces 
pressure  transients.  Ir.  this  regard,  stress  waves  can  trigger  rapid  gasifi¬ 
cation  of  the  melt  layer  by  partially  collapsing  cracks.  Such  stress  waves 
originate  w.ier.  cracks  propogate  into  relatively  high-pressure  cavities.  Sub¬ 
sequently,  their  amplitude  may  be  increased  by  reflection  of  the  stress  wave 
by  stiff  media  such  as  a  motor  case. 

The  effect  of  stress  waves  is  two  fold.  The  first  is  to  raise  the  gas 
pressure  by  reducing  the  void  volume  and  thereby  increase  the  rate  of  heat 
transfer  into  the  burning  in  an  propellant..  In  turn  the  increased  heat  fluxes 
accelerate  the  burning  in  an  exponential  fashion.  The  second  effect  of  stress 
waves  is  to  restrain  crack  expansion  while  the  pressures  always  peak  before 
decaying.  When  sufficiently  strong  stress  waves  act  upon  burning-propellant 
cracks  with  adequate  melt  (depends  upon  stress  wave),  pressure  rises  of  the 
oeder  of  tens  of  kbars  can  be  generated  within  a  fraction  of  a  usee  (See 
Section  4).  This  statement  is  conservative  in  that  it  neglects  ejection  of 
melt  and/or  particulates  into  the  hot  combusion  gases  or  other  processes  that 
apparently  speed  consumption  of  the  melt  (See  Section  5.3).  Step-wise  pres¬ 
sure  waves  with  amplitudes  of  the  order  of  tens  of  Kbars  are  known  to  initiate 
a  composite  propellant  containing  HMX  in  impac.e  experiments.6 

One  can  conceive  of  three  situations  in  which  adequate  melt  layers  needed 
for  DDT  can  develop.  These  are: 

1.  burning  in  cracks  or  debonds  at  near-ambient 
pressures  wherein  melt  layers  ;re  greatest 
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2. 


sustained  low  rates  of  heating  (less  than  that 
associated  with  burning  or  1  bar  of  pressure)  of 
unignited  propellant  surfaces, 

3.  accumulation  of  flowing  melt  within  cracks 

Of  the  three  only  the  third  appears  placeable.  In  this  regard,  Item  1  is  highly 
unlikely  in  that  it  requires  burning  under  conditions  of  sustained  near-ambient 
pressures.  Moreover,  Item  2  appears  unlikely  in  that  it  is  difficult  to  com¬ 
prehend  how  low  rates  of  heating  could  be  sustained  long  enough  to  develop 
substantial  melt. 

Item  3  seems  most  likely  because  of  the  appreciable  velocities  of  qases 
flowing  into  cracks  from  high-pressure  cavities  (i.e.  of  the  order  of  hundreds 
of  meters/sec),  and  the  relatively  slow  rate  of  gasification  of  newly  melted 
HMX  (i.e.  more  than  an  order  of  magnitude  slower  than  molten  HMX  during  burn¬ 
ing).  Melt  will  accumulate  if  sufficient  molten  HMX  is  moved  over  the  crack 
surfaces  by  the  gas  flows  rather  remaining  stagrent  or  being  swept  into  the 
gas  stream.  Questions  regarding  the  latter  remain  to  be  resolved.  If  melt 
flows  over  th  ->ck  surfaces,  then  the  cooler  most  recent  melt  can  mix  with 
hotter  melt  e  crack.  The  result  would  be  a  lowering  of  the  Tatters  tem¬ 
perature  and  nence  its  rate  of  gasification.  By  this  process,  substantial 
melt  may  accumulate  even  though  the  propellant  Is  burning  at  high  pressures 
wherein  the  melt  layer  Is  normally  very  thin.  Only  by  this  process,  does  it 
appear  possible  to  generate  the  amounts  of  melt  needed  for  DDT  when  a  single 
crack  or  debond  Is  present. 

Appreciably  less  melt  is  required  for  DDT  when  several  burning  propellant 
cracks  are  present.  The  latter  is  because  the  consequence  of  stress  waves 
acting  upon  a  burning-propellant  crack  are  stress  waves  of  greater  amplitude. 

By  the  process  progressively  stronger  stress  waves  can  be  generated  as  they 
move  from  one  burning  crack  to  other  similarly  oriented  burning  cracks.  Tran¬ 
sient  gas  pressures  will  progressively  increase  In  successive  cracks  and  read' 


their  peaks  in  shorter  and  shorter  times.  In  this  regard,  rapid  pressure-rls 
times  are  an  important  requirement  for  DDT  as  well  as  high-amplitude  pres  es 
of  the  order  cited  earlier.  That  is  because  the  pressure  waves  must  dev 
extremely  steep  formats  to  cause  detonation.  At  pressures  in  excess  of  -jhly 
10  kbar,  higher  pressure  portions  of  the  wave  will  eventually  catch  •  .i 
the  lower  pressure  portions  of  the  wave  provided  sufficient,  pro  _  pre¬ 

sent.  That  Is  why  rapid  pressure  rises  are  important. 
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6.2  EXPERIMENTS 


6.2.1  MeZt  Laytx 

The  mass  of  the  melt  or  foam  present  during  steady  burning  of  the  HtiX- 
based  composite  propellant  was  determined  at  ambient  pressure.  Its  purpose 
was  two  fold;  namely  to  determine  whether  or  not  there  is  sufficient  melt  to 
produce  ODT  experimentally  and  to  gain  a  better  appreciation  of  the  validity 
of  the  substantial  melt  predicted  for  pure  HMX  at  1  bar  (See  figure  4). 

The  determinations  involved  calorimetric  measurements  of  the  total  heat 
content  in  burning  propellant  specimens,  and  thermocouple  measurements  of  the 
heat  content  in  the  solid  propellant.  From  these  measurements  the  heat  content 
of  the  melt  was  obtained  and  divided  by  the  enthalpy  of  pure  HMX  at  its  reaction 
temperature  to  estimate  the  mass  of  the  melt  layer.  It  was  found  to  be  in 
substantial  agreement  with  that  predicted  for  pure  HMX.  Computer  runs  indica¬ 
ted  it  was  adequate  to  produce  ODT  by  the  experimental  design  described  in 
Section  5  of  this  report. 

6.2.2  WT  ExpeJusmnti, 

A  total  of  eighteen  experiments  were  conducted  in  an  attempt  to  initiate 
a  burning  HMX  based  propellant.  They  involved  partial  closure  of  the  void 
space  above  burning  propellant  in  a  manner  analogous  to  that  produced  by  stress 
waves  acting  upon  burning  propellant  cracks.  The  latter  was  achieved  by  driv¬ 
ing  a  piston  into  a  chamber  containing  a  burning  cylinder  of  the  propellant. 

In  each  experiment,  partial  closure  of  the  void  space  over  the  burning 
propellant  cylinder  resulted  in  blowing  the  piston  and  driver  weights  out  of 
the  chamter.  Detonation  did  not  occur.  Lack  of  detonation  was  attributed  to 
inadequate  closure  of  the  void  space  predicted  for  DDT.  There  are  three  pos¬ 
sible  causes  of  t.ie  lack  of  adequate  closure.  These  are: 

•  excessive  drag  on  piston  caused  by  contaminants 
being  weighed  between  piston  and  chamber  walls 

•  thermal  deformation  of  chamber  walls  exposure  to 
combustion  gases 

•  more  rapid  consumption  of  melt  or  foam  layer  than 
anticipated  theoretically 
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Several  materials  were  used  to  lubricate  the  bearing  surfaces.  Initial 
lubricants  were  not  satisfactory  In  that  the  piston  was  slowed  by  passage  of  com¬ 
bustion  products  into  the  0.001  clearance  between  the  piston  and  chamber.  Sub¬ 
sequently,  closure  was  improved  using  a  film  of  silicon  grease  of  thickness 
less  than  0.001  in.  over  which  produced  graphite  was  sprinkled.  The  possibility 
of  piston  stoppage  due  to  thermal  deformation  of  the  chamber  walls  was 
ascessed  by  increasing  the  piston/chamber  clearance  several  fold.  Lack  of 
closure  Improvement  suggested  thermal  distortion  was  the  cause  of  the  problems. 
From  these  results,  it  appears  that  Inadequate  closure  was  due  to  extremely 
rapid  consumption  of  the  melt  or  foam  layer.  This  conclusion  is  consistent 
with  a  two-fold  reduction  of  the  void  space  when  a  vent  hole  was  pro¬ 
vided  to  lessen  the  pressure  rise.  Two  possibilities  may  be  advanced  of 
for  extremely  rapid  burning.  The  first  is  blow  off  of  melt  and  particulate 
into  the  hot  gas  stream  wherein  they  are  most  rapidly  gasified.  The  second 
is  Inadequacy  of  the  expression  pQ^  used  to  described  the  internal  heating. 

This  expression  implies  a  constant  quantity  of  the  heat  is  transferred  to  the 
foam  layer  per  unit  mass  of  gas  evolved.  Residence  times  of  the  gases  in  the 
foam  will  affect  the  extent  to  which  the  gases  decompose  In  the  foam  and  hence 
the  amount  of  heat  transferred  to  the  foam.  In  highly  dynamic  burning  under 
consideration,  order  of  magnitude  variations  of  residence  time  may  be 
expected---fi rst  due  to  order  of  magnitude  changes  in  the  amount  of  foam 
present,  and  second  due  to  reductions  of  gas  bubbles  brought  about  by  rising 
pressures. 

The  above  conclusions  pertaining  to  consumption  of  the  melt  layer  do  not 
alter  the  conclusion  presented  in  the  report  in  regard  to  DDT.  Instead  they 
make  DDT  More  likely  than  anticipated  by  speeding  gasification  of  the  foam 
provided  adequate  crack  close  is  achieved. 

6.2.3  Rec  ormendatcotu 

Based  upon  conclusions  drawn  from  Section  6.2.2,  appreciably  higher  pis¬ 
ton  speeds  and/or  smaller  void-volume  heights  are  needed  to  validate  our  DDT 
theory.  It  is  anticipated  that  piston  speeds  at  least  several  times  greater 
than  those  used  herein  should  be  used  depending  upon  the  mass  of  the  piston 
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and  Its  appendages.  Driving  the  piston  with  gas  pressure  would  provide  a 
controlled  versatile  means  for  developing  the  necessary  momentum.  In  addi¬ 
tion  It  Is  desirable  to  reduce  the  Intlal  void  space. 

The  above  device  Is  also  ameanable  to  experiments  with  which  to  assess 
how  substantial  melt  layers  are  consumed  during  accelerated  burning.  Here 
pressure  transducers  would  be  needed  along  with  means  to  determine  the  void 
space  dynamically.  A  comparison  of  the  dynamic  pressures  with  analytic  pre¬ 
dictions  would  serve  to  Indicate  the  adequacy  of  existing  burn  theories  and 
means  for  improvement  thus  predictions  capability. 

While  the  above  requirements  are  rather  costly,  it  Is  believed  that  they 
should  shed  considerable  light  on  DDT.  Because  of  the  extremely  dynamic  burn 
conditions,  they  would  also  severely  test  existing  theories  of  how  such  pro¬ 
pellants  burn  dynamically. 
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APPENDIX  A 


METHOD  FOR  PREDICTING  TRANSIENT  TEMPERATURE  IN 
SOLID-PHASE  OF  PROPELLANTS 

In  order  to  predict  dynamic  burning  of  propellants,  it  is  necessary  to 
compute  the  propellant  temperatures  under  rapidly  changing  pressures  and 
conditions.  In  this  regard,  finite  difference  methods  requires  periodic 
refinment  of  various  spacial  increments  when  the  burning  becomes  extremely 
dynamic.  As  a  consequence  such  methods  are  not  computationally  efficient  in 
dealing  with  problems  in  which  the  temperature  gradients  within  the  propellant 
become  extremely  prounounced. 

In  view  of  the  above,  a  method  of  sources  and  sinks  was  developed.  It 
predicts  transient  temperatures  within  solid  propellants  at  desired  depths  and 
times.  The  method  may  be  applied  to  burning  propellants  with  or  without  a 
melt  or  foam  layer  on  the  propellant's  surface. 8 

1.  DESCRIPTION  OF  METHOD 

In  this  method,  the  propellant  is  assumed  to  be  semi-infinite  in  that 
heat  penetrates  only  relatively  shallow  depths  during  burning.  The  mputa- 
tional  method  utilizes  a  fixed  cartesian  coordinate  system  in  which  i.,e  rnelt  or 
burn  interface  is  considered  to  move  in  a  step-wise  fashion  with  respect  to 
time  as  shown  in  Figure  A-l.  As  will  be  shown  later,  mass  removal  is  accounted 
for  by  adjustment  of  the  heat  fluxes. 

Prior  to  the  start  of  melting  or  burning  the  "moving  boundary"  is  of 
course  stationary.  During  subsequent  time  steps  ,*t^,  the  interface  is  con¬ 
sidered  to  move  from  depth  x^_,  to  x^ .  Time  average  heat  fluxes  crossing  the 
moving  solid-propellant  boundary  during  &t^  are  represented  by  q...  These 
fluxes  are  evaluated  from  the  boundary  conditions  at  the  moving  boundary,  and 
vary  with  the  propellant,  and  pressures  and  temperatures  to  which  it  is 
subjected.  Each  flux  q.  or  more  precisely  is  located  at  particular  depth 
that  yields  the  same  quantity  of  heat  entering  the  propellant  at  x^  during 
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Oesired  Flux 


Boundaries  Prior  to  Melting 
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Figure  A-l.  Schematic  of  computational  procedure. 
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Atj  as  would  be  produced  by  the  moving  boundary.  Means  for  determining  the 
appropriate  depth  are  described  in  Section  1.2  of  this  appendix. 

In  this  discussion  to  follow,  the  subscript  i  designates  values  of  q,  t 
and  x  associated  with  prior  time  steps  while  the  subscript  j  Indicates  values 
for  the  most  recent  time  step  At-.  First  we  shall  discuss  means  for  achieving 

J 

the  given  flux  q.  at  the  desired  depth  x.’. 

J  * 

1.1  APPLIED  HEAT  FLUXES 

The  consequence  of  applying  fluxes  at  a  depth  x^'  during  Ati  is  to  create 
conductive  fluxes  q.  •  at  the  depth  x^'  during  the  most  recent  time  step  At.. 

■  f  J  J  J 

These  conductive  fluxes  must  be  eliminated  to  account  for  propellant  removal. 
Therefore  in  order  to  produce  a  desired  flux  q.  at  x.',  it  is  necessary 

J  J 

to  apply  a  flux  q.'  equal  to  the  difference  between  q.,  and  the  undesirable 

J  J 

conductive  flux.  To  illustrate  the  above  mathematically  we  shall  represent 
the  undesirable  time-average  conductive  flux  produced  at  x-'  during  At-  by 

J  J 

each  of  the  prior  fluxes  q . ‘  and  q. ..  Thus,  the  time-average  flux  q.'  to  be 

•  I  J  s) 

applied  at  x s *  during  At.  is 

J  J 


j  -  1 


i  =  1 


Clearly  one  must  know  q.j  in  order  to  determine  the  flux  q ^ 1  that  needs 
to  be  applied  at  x-  to  achieve  the  desired  flux  q..  For  this  purpose,  consider 

s)  V 

the  temperature  rises  produced  within  a  semi-infinite  body  of  uniform  initial 
temperature  by  flux  q,'  of  constant  magnitude  applied  at  a  depth  x-'  for  all 

J  J 

times  to  starting  at  t-  ,.  The  temperature  rises  beneath  x.'  are  given  by:9 

J"*  J 


(A-2) 


where  C  =  2//KpUj^  and  t  _>  tj_;  =  At^  +  At2  •  .  •  At-_j.  Equation  A-2  applies 
to  the  most  recent  flux  ' .  In  order  to  determine  the  temperature  rises 
produced  by  prior  fluxes  q^ one  must  recognize  that  temperature  rises 
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produced  by  a  series  of  constant  fluxes  are  additive.  Replacing  the  subscript 
j  In  the  above  expression  by  i  yields  the  temperature  rise  produced  by  the 
flux  q^  *  applied  for  all  times  t  greater  than  t-^.  Replacing  t^  by  t^ 

yields  the  temperature  rise  produced  by  the  same  flux  when  It  is  applied  for 
all  times  t  greater  than  t..  Subtracting  the  latter  temperature  rises  from 
the  form  yields  the  temperature  rises  produced  by  q. '  applied  over  the  time 
step  from  t;  to  t. .  The  result  is  given  below: 

1-1  I 


c  V 


x  -  xi 


/t  -  t4  i erfc  - 

1-1  2  r^TCt  t.  ) 

i  ■  i 


-  /t  -  t.  ierfc 


x  -  x. 


2  Mi  -  tT 


(A-3) 


Here  x  >  x^ 1  and  t  >  t^ . 

Multiplying  the  above  expression  by  the  negative  value  of  the  propellant's 
thermal  conductivity  K  and  differentiating  with  respect  to  x  yields  the  time- 
dependent  fluxes  at  x  presented  below 


x  -  x. 


x  -  x. 


erfc 


2  /a(t  -  tj_  ) 


-  erfc 


2  /o(t  - 


(A- 4) 


The  subscript  i  ranges  from  1  to  j  -  1. 

To  find  the  time-average  fluxes  q, ■  at  x^  over  the  time  step  At.:,  we 

1  J  J 

shall  replace  x  and  t  of  Equation  A-4  by  x.-'  and  t-  +  t,  respectively.  The 

sj  s)  1 

integral  of  the  result  over  the  time  step  divided  by  Atj  yields  the  desired 
time-average  flux  q^j  presented  below: 


■u 


At 


j 


erfc 


1L 


-  xi 


erfc 


-  X, 


2  /oTT:  -  ti  +  t) 
j*  i  i-i 


2  M^~i  -  t.  +tT 

J-  1  ! 


dr 
(A-5) 


As  one  would  expect  q^j  is  zero  if  the  propellant  interface  does  not  recced, 
i.e.,  x.'  =  xi '  =0.  Performing  the  integration  of  Equation  A-5  yields; 

J  * 
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qu  =  )  -  F(y *)  -  +  Hy>)] 

3 


(A-6) 


Where  the  function  F  is  expressed  below  in  terms  of  the  normal  distribution 
function  represented  by  N. 


F(y)  =  -  +  y2  ^  exp  (-y2)  +  4  Nf^Zy) 


and  y! ,  y2,  y3,  and  y4  are 


y,  •  Ixj'  -  x.')/2 


y2  =  (Xy1  -  x. '  )/2  /a(t,  '  t,  no 

J  •  J  i  i  - 1  j 

y3  *  (Xy'  -  X.  *)/2  -  t.  ) 

J  I  J-l  1 


y4  =  (Xj'  -  x^/Z  /a(tj_Y  -  t.  +  AtjJ 
Substituting  q . ,  into  Equation  A-l  yields  the  flux  q.‘  that  must  be 

'  J  %) 

applied  at  x^'  during  Atj  in  order  to  achieve  the  net  desired  flux  q^.  Means 
for  defining  the  depth  x,‘  are  described  in  the  following  section. 

J 


1.2  DEPTHS  AT  WHICH  FLUXES  ARE  APPLIED 


i  i 


I  > 

1  i 


The  depth  x,'  at  which  q^'  is  applied  affects  the  heat  transfer  into  the 
propellant  at  the  depth  xA  Location  of  the  depth  x:'  becomes  more  critical 

w  1 J 

as  the  size  of  the  spacial  increment  Ax.  Increases.  For  optimal  utility,  it 

V 

Is  desirable  that  the  method  be  capable  of  using  relatively  large  Ax^  wherein 

J 

the  depth  x., '  Is  Important.  The  desired  depth  x,'  may  be  defined  in  terms  of 

J  J 

the  fraction  £  of  the  displacement  Axj  during  At^  as  follows: 


V  =  xj-.  *  s  (xj  •  xj-,>  *  xj-,  +  5  ixj 


IIT  RESEARCH  INSTITUTE 


(  A~7  ) 


70 


T  '  r  - 


f  ’  •* 


] 


The  factor  £  varies  with  the  acceleration  of  the  interface  and  the 
thermal  conduction  parameter  3  described  below: 

i 

e  *  AXj/(2  (A-8) 


Two  sets  of  calcualtions  were  used  to  determine  the  above  dependence. 

First  a  given  displacement  Axj  was  subdivided  into  small  increments 
wherein  the  location  of  the  applied  times  ceased  to  be  important.  Then  the 
quantity  of  heat  passing  the  depth  x,  during  At.  was  computed  by  su.nning  the 

w  w 

time-dependent  fluxes  predicted  by  Equation  A-4.  It  is  this  quantity  of  heat 
that  should  be  effected  by  the  correct  value  using  the  single  spacial 
increment  Ax-. 

s) 

The  final  set  of  calculations  involved  varying  £  with  the  single  incre¬ 
ment  Ax-  until  the  heat  entering  the  depth  x.  during  At,  equaled  that  predicted 

J  J  J 

above  using  extremely  fine  spacial  and  temporal  increments.  Results  are 
presented  in  Table  A-l. 


TABLE  A-l.  £  VALUES 


Changes  of  Rate  of  Melting 

Ouring  Time  Step  At,,  £  Values  (dimensionless) 

percent  J  3  =  0/(54  6  =  0.  OS  6  =  07l  2 


-80 

0.289 

0.295 

0.300 

-60 

0.321 

0.325 

0.329 

-40 

0.322 

0.335 

0.339 

-20 

0.336 

0.339 

0.341 

0 

0.335 

0.337 

0.341 

40 

0.334 

0.336 

0.338 

80 

0.329 

0.331 

0.332 

150 

0.319 

0.320 

0.321 

300 

0.297 

0.299 

0.300 

600 

0.270 

0.270 

0.270 

1.3  TEMPERATURE/FLUX  PREDICTIONS 

Temperature  rises  produced  by  the  heat  fluxes  q ^ '  and  q^ '  at  depths  x^' 
and  Xj '  over  time  steps  Atj  and  At^  are  given  by  Equations  A-2  and  A-3 
respectively.  Values  for  q  ^ '  or  '  are  determined  from  Equations  A-l  and 

A-5  while  the  depths  x^'  and  x,. '  are  found  as  described  by  the  previous  section 
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*1'  i  7r'* 


In  that  the  temperature  contributions  are  additive,  the  temperature  of  the 
surface  of  the  solid  propellant  at  the  melt  or  burning  interface  at  the  end  of 
the  time  step  Ati  is  given  by  Equation  A- 8. 

J 


j  -  1 

T<Xj)  ■  C  £  y 
i  =  1 


VA 


v't.  -  t,.  ierfc 

2  ^Tt"“t~T 


xj  '  xi 


A;  -  t .  ierfc  - 

J  1  2  Mi —  tjl 


(A-8) 


X  •  —  X  • 

+  C/At~  q  ’  ierfc  — - i  +  T 

J  J  2  w^AtJ  0 


For  propellants  that  melt,  T(Xj)  equals  the  melt  temperature  Tm  once  melting 
starts. 

The  heat  flux  entering  the  solid  propellant  at  x.,  t.,  namely 

J  J 


J 

i 


-  1 


I 


=  1 


xi  ' 

erfc  - J . . . 

.  2  '"ltj  •  v,1 


-  erfc 


xj  *  xi 


2  Mi i 


( A-9) 


xi  -  x. 

q,*  erfc  -  ^ - - 

2  /a(t.  - 


equals  the  flux  supplied  externally  when  the  propellant  is  not  melting.  When 
propellants  melt,  it  is  necessary  to  combine  the  fluxes  heat  flux  expressed 
above  with  expressions  describing  the  external  and  internal  generated  heating 
of  the  melt  layer  to  arrive  at  the  temperature  of  the  melt  and  its  rate  of 
gasification.  The  latter  expressions  are  presented  in  Section  2. 


2.  VALIDATION  OF  METHOD 


2.1  PROPELLANT  WITH  FOAM  LAYER 

Here  we  are  concerned  with  checking  the  method  by  assessing  how  well 
dynamic  burning  approaches  known  steady  burning  under  fixed  external  conditions. 
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Initially,  the  propellant  is  considered  at  a  uniform  ambient  temperature. 
Thereafter  the  propellant  Is  exposed  to  a  constant  heat  flux  equal  to  that 
present  during  steady  burning  associated  with  a  preselected  constant  pressure. 
Dynamic  calculations  were  then  conducted  until  the  changes  of  the  bum 
velocities,  foam  temperatures,  and  foam  masses  were  insignificant.  These 
results  are  presented  in  Table  A-2  under  the  columns  titled  Model  Predictions. 
Values  predicted  by  Equations  13  and  14  are  in  the  colunms  titled  Analytical 
Predictions. 

It  may  be  observed  that  the  model  predictions  are  in  good  agreement  with 
the  analytical  predictions.  They  indicate  that  the  cunmulation  of  errors  by 
the  numerical  calculations  is  relatively  small. 

?ficpdtZay\A  (tlcthotU  Foam  LaytA 

In  this  case  validation  was  achieved  by  comparing  the  regression  rate 
predictions  with  those  obtained  by  Kooker  using  finite  differences.  The 
calculations  involved  the  KTSS  combustion  model.4  Initially,  Kooker10  con¬ 
sidered  the  solid  propellant  burns  at  a  steady  regression  rate.  Thereafter 
the  pressure  V  (dimensionless)  is  considered  to  increase  as  follows 

P  =  1  +  2.5  (1  -  exp  (-2.5  t))  (A-10) 

where  T  represents  the  dimensionless  time. 

The  consequences  of  increased  pressure  is  increased  regression  rates. 
Figures  A-2  and  A-3  present  the  results  of  Kooker10  following  the  increase  of 
pressure  described  by  Equation  A-10.  In  these  figures,  N  and  M  are  pressure 
and  temperature  rise  exponents  describing  steady  and  nonsteady  burning  rates, 
respectively.  H  is  proportional  to  the  internal  heating.  Our  results  are 
presented  by  dots. 

It  may  be  observed  that  our  regression  rate  predictions  are  in  reasonably 
good  agreement  with  those  obtained  by  Kooker10  considering  both  methods  are 
numerical.  Number  of  time  steps  used  to  generate  the  IITRI  results  of  Figures 
A-2  and  A-3  were  73  and  188,  respectively,  and  involved  execution  times  of 
5  and  37  seconds  on  a  1109  computer. 
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Figure  A- 2. 


Vatidation  of  TITRI  method  H  -  0.75 
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APPENDIX  6 


DETERMINATION  OF  CONSTANTS  d  AND  c2  FOR 
HEAT-TRANSFER  COEFFICIENT  h 

Here  we  shall  determine  the  constants  ci  and  c2  of  Equation  7  by  use  of 
estimated  foam  masses  during  steady  burning  of  HMX.  To  this  end,  two 
equations  must  be  solved  for  Ci  and  c2. 

The  first  equation  Is  based  upon  Equations  6  and  7  for  steady  burning 
as  follows: 

qp  =  cl(exp(-E/Tf(P)))C2(Tf  -  TJ  (6-1) 

To  support  steady  burning,  qp  must  also  satisfy 

\  •  »?ftCp(Tm  -  To)  f  Qm3,  (B-3) 

Equating  the  above  expressions  for  qp  and  using  Equation  13  yields 

c,(exp(-E/Tf(P)))C2(Tf  -  TJ  =  paPn(Cp(Tm  -  Tq)  +  QM)  (B-3) 

Equation  B-3  represents  the  first  of  the  two  equations. 

The  secord  equation  Is  obtained  from  Equations  13  and  14  for  steady 
velocities  r^.  It  is  given  by 

paPn  «  Mf(P)  Z  exp(-E/Tf(P) ))  (B-4) 

At  any  given  pressure  P,  all  of  the  parameters  of  Equations  B-3  and 
B-4  are  known  except  for  the  constants  ci  and  c2,  and  the  foam  mass  and 
temperature  at  the  given  pressure.  To  determine  the  constants  Ci  and  c2 
it  is  necessary  to  know  the  melt  mass  M^(P)  at  two  pressures  P.  In  this 
regard,  Boggs  has  photographed  the  "frozen"  HMX  foam  layer  following  ex¬ 

tinguishment  of  the  burning  by  rapid  pressure  relief.  Bogg's  photographs 

Indicate  that  the  steady-state  foam  thickness  at  34  bars  is  roughly  2D  urn 

IIT  RESEARCH  INSTITUTE 


77 


thick;  at  68  bars  It  appears  to  bo  about  half  as  thick  as  that  at  34  bars  of 
pressure. 

Assuming  that  the  density  of  the  "frozen  foam"  Is  half  the  density 
(1.9  g/cms)  of  solid  high-density  HMX  propellant  yields 

Mf{34)  =  0.0024  g/cm2  (B-5) 

Rf ( 68 )  =  0.0012  g/cm2  (B-6) 

At  this  point,  one  may  raise  a  question  regarding  the  correspondence 
between  the  mass  of  "frozen"  foam  with  the  foam  mass  present  during  steady 
burning.  In  other  words,  how  much  does  the  foam  mass  change  following 
pressure  relief  due  to  further  outgasslng  and  melting? 

The  one-dimensional  model  (foam)  was  used  to  resolve  the  above  question. 

It  Indicated  that  the  foam  mass  at  34  bars  will  decrease  by  13.4  percent 
during  pressure  relief  while  the  foam  mass  at  68  bars  will  decrease  by  19.4 
percent.  The  somewhat  greater  reduction  (percentagewise)  of  the  foam  mass 
present  at  68  bars  Is  due  to  Its  higher  temperature.  The  result  Is  more 
rapid  mass  loss  of  Initially  smaller  amounts  of  foam. 

While  the  above  mass  losses  are  significant,  they  are  probably  small 
compared  to  errors  in  estimating  the  thickness  and  density  of  the  "frozen" 
foam.  For  this  reason,  the  foam  masses  presented  by  Figure  4  shall  not  be 
corrected  until  better  measurements  of  the  "frozen"  foam  mass  are  available. 
Substituting  each  pair  of  M^,  P  values  of  Equations  B-5  and  B-6  Into  Equation  B-4 
and  solving  for  T^(P)  with  the  propellant  properties  given  in  Appendix  F 
yields 


Tf  (34 )  =  689. 5°K 
Tf (68)  =  713. 0°K 


(B-7 ) 
(B-8) 


Substituting  the  above  T  p 

e,Mt,0ns  '"•'"H  the  two  Jk'nowns  c'Sa„dnf°jEqUat,°"  B-3yU'“  *» 


1  cal/cm2-sec-°K 

c2  *  0.338  (dimensionless) 


(B- 9) 
(8-10) 
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APPENDIX  C 
IGNITION  CRITERIA 

Ignition  refers  to  the  start  of  propellant  burning.  Knowledge  of  when 
Ignition  occurs  Is  Important  In  that  the  heat  flux  undergoes  rapid  change 
following  Ignition. 

Ignition  represents  a  complex  phenomenon  resulting  In  exponential 
temperature  rises  of  the  evolved  reaction  gases  caused  by  exothermic  decom¬ 
position.  Ignition  Is  a  function  of  a  number  of  factors  foremost  of  which  are: 

•  temperature,  composition  and  flow  rates  of  the 
gases  evolved  by  the  heated  propellant,  and 

•  physical  and  thermal  environment  into  which 
the  gaseous  products  are  discharged. 

Fortunately,  as  will  be  seen  later  in  this  appendix,  one  need  not  conduct 
detailed  analyses  of  each  of  the  above  phenomena  to  approximate  Ignition 
times.  Nevertheless  one  should  be  aware  of  the  phenomena  In  designing  ex¬ 
periments  with  which  to  establish  Ignition  criteria  for  the  propellant/ 
conditions  of  Interest. 

At  present.  Ignition  criteria  are  based  upon  propellant  temperature,  or 
combinations  of  a  constant  incident  flux  and  duration.12  The  former 
criterion  Is  of  greater  value  in  that  it  applies  to  time- dependent  heat 
fluxes.  Ignition  Is  predicted  wnen  the  temperature  of  the  propellants  surface 
or  at  some  propellant  depth  exceeds  a  critical  temperature. 

In  the  remainder  of  this  section  we  shall  examine  the  effect  of  dis¬ 
regarding  gas  evolution  upon  the  Ignition  time.  Before  doing  so  It  should 
be  noted  that  a  one-to-one  correspondence  does  not  exist  between  foam  (or 
surface)  temperature  and  rates  of  gas  evolution. 

Table  C-I  presents  predicted  foam  temperatures  following  exposure  of 
HMX  to  specified  Incident  heat  fluxes.  Internal  heating  Is  also  provided  for 
even  though  It  is  not  specified  in  the  table.  The  Initial  temperature  of  the 
propellant  Is  294°K.  Foam  or  propellant  surface  temperatures  are  specified 
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lit. _ 


at  times  at  which  the  gas  evolution  rates  equal  specific  fractions  6  of  the 
steay  burning  rate  of  0.03  cm/sec  at  1  atomsphere  of  pressure. 

TABLE  C-l.  PREDICTED  FOAM  TEMPERATURES  AND  REGRESSION  RATES  OF  HEATED  HMX 


Temperature  T^  yielding  to  specified  fraction  6 

Incident  Heat _ _ ofTthe  steady  burning  rate,  KJ _ 

FI  ux* 


cal/cm2-sec 

6  =  0 

6  *  0.1 

6  =  0.5 

6  =  0.9 

6  =  1.0 

0.5 

555.0 

565.6 

580.9 

587.7 

588.9 

(211.44) 

(237.19) 

(239.33) 

(239.58) 

(239.60) 

0.1 

555.0 

570.30 

585.9 

589.1 

592.3 

(52.86) 

(61.99) 

(64.45) 

(64.71) 

(64.74) 

2.0 

555.0 

579.0 

588.9 

594.9 

596.1 

(13.25) 

(16.50) 

(18.25) 

(18.50) 

(18.53) 

10.0 

555.0 

603.5 

607.8 

611.6 

612.4 

(0.529) 

(0.668) 

(0.875) 

(0.959) 

(0.974) 

30.0 

555.0 

629.1 

630.4 

631.8 

632.2 

(0.0587) 

(0.0697) 

(0.0902) 

(0.01031) 

(0.1057) 

100.0 

555.0 

665.6 

664.5 

664.5 

664.6 

(0.00529) 

(0.00577) 

(0.00681) 

(0.00760) 

(0. 00770) 

♦Incident  heat 

fluxes  are  from 

some  external 

heat  source 

Times  of  occurrence  given  In 

parentheses  in 

seconds . 

First 

It  should  be  observed  that  more 

!  time  is  needed 

to  initiate  melting 

than  the  remaining  time  needed  to  achieve  steady-state  rates  of  gasification 
(5  =  1.0).  The  latter  is  due  to  rapid  decompostion  of  HMX  at  and  above  Its 
high  melt  temperature. 


From  Table  C-l  It  may  be  observed  that  an  ignition  crite.'ion  based  solely 
upon  a  glvan  temperature  implies  that  Ignition  occurs  with  differing  rates  of 
gas  release.  The  latter  reemphasizes  the  question  raised  earlier  regarding  the 
significance  of  the  differing  rates  of  gas  evolution  upon  the  Ignition  time. 
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In  order  to  perform  the  above  assessment  we  shall  assume  that  the  i 

critical  Ignition  temperature  lies  between  the  melt  temperature  of  555  K  and 
the  lowest  temperature  cited  in  Table  C-I  for  6=1.  A  6  value  of  1  Is  chosen 
since  It  corresponds  to  a  regression  rate  equal  to  the  steady  burning  rate. 

Based  on  the  above  assumptions,  the  critical  temperature  lies  between  555.0 
and  588.9  K. 

Applying  the  above  temperature  range  to  Table  C-l  indicates  that  the 
Ignition  time  Is  between  211.4  and  239.6  sec  when  the  Incident  flux  Is  0.5 
cal/cm2-sec;  and  between  0.00559  an  a  linearly  Interpolated  value  of  0.00544  j 

sec  when  the  incident  flux  is  100  cal/cm2-sec.  Notice  that  tie  above  times  i 

differ  by  only  13.3  and  2.8  percent  for  the  two  fiuxes.  These  results  sug-  j 

gest  that  ignition  temperature  provides  a  reasonably  accurate  means  for  j 

predicting  Ignition  at  least  for  HMX.  Rates  of  gas  evolution  are  of  secondary  j 

Importance.  ] 

For  the  present  we  shall  assume  that  the  ignition  temperature  of  HMX  j 

1 

equals  the  mean  temperature  (572°K)  associated  with  the  range  of  temperatures  j 

cited  earlier.  Ignition  times  t^  required  to  achieve  a  foam  temperature  of  j 

572°K  may  be  determined  from  Table  F-l  as  a  function  of  q.  They  are  given  j 

approximately  by  ; 

j 

t.  =  59/q2  (C-l)  ; 

where  t^  Is  In  seconds  when  q  is  in  cal/cm2-sec.  Experimental  results  pre¬ 
sented  In  reference  12  Indicate  that  the  exponent  2  of  Equation  F-l  can  vary  j 

from  about  1.6  to  2.0  depending  upon  the  propellant.  Most  propellants  Involve  j 

coefflcents  only  a  fraction  of  the  coefficient  59  presented  by  Equation  F-l. 

The  relatively  large  coefficient  is  attributed  to  the  high  melt  temperature  of 
HMX  and  the  fact  that  much  of  the  "Ignition  time"  is  expended  In  Initiating  j 

melting.  In  this  regard,  times  to  Initiate  melting  are  proportional  to  the  ; 

square  of  the  difference  between  the  melt  temperature  and  the  initial  tem¬ 
perature  of  the  propellant.9 
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APPENDIX  D 


DYNAMIC  REGRESSION  RATES  OF  PROPELLANTS 
HAVING  NO  FOAM  LAYER 

Dynamic  regression  rates  of  completely  solid  propellants  are  ususally 
predicted  in  terms  of  the  surface  temperature  of  the  propellant.3  No  account 
is  made  for  propellant  temperatures  beneath  the  surface  or  of  the  mass  of 
propellant  involved  at  elevated  temperatures.  Clearly  the  greater  the  amount 
of  propellant  at  elevated  temperatures,  the  greater  the  rate  of  gasification 
should  be.  It  is  with  this  reason  that  the  following  analysis  is  conducted. 

To  account  for  variations  in  the  propellant  temperature  as  a  function  of 
depth,  dynamic  regression  rates  are  predicted  by  integrating  the  Arrenhuis 
relationship  as  follows: 


r 


f  * 


exp  (-E/(T(x)  +  T0))  dx 


(D-l) 


where  x  represents  depth  beneath  the  bum  surface  at  time  t.  Due  to  the 
presence  of  pronounced  temperature  gradients  during  burning.  Equation  D-l  need 
only  be  Integrated  over  shallow  depths. 

A  simpler  more  approximate  expression  for  r  may  be  achieved  by  approxi¬ 
mating  the  temperature  rises  by 

T(x)  =  T(o)  exp  -  (C,x)  (D-2) 

where  C3  is  given  by 

C3  =  q/ (  KT( o ) )  (D-3) 

Substituting  T(x)  of  Equation  D-2  into  Equation  D-l,  and  integrating  by 
parts  with  the  neglect  of  second-order  terms  yields 
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(D-4) 


rf  a,  ^  exp  (-F/(T(o)  +  T0))  (T(o)  +  TQ)2 

q£ 

It  may  be  observed  that  the  above  expression  differs  in  form  from  that  con¬ 
ventionally  used.3  Usually  the  exponential  function  is  multiplied  by  a  con¬ 
stant.  Here  it  includes  the  variable  factor  (T(o)  +  TQ)J/q  to  account  for 
variations  in  the  depths  of  heated  propellant. 

In  order  to  use  Equation  D-4,  it  is  necessary  to  determine  the  time- 
dependent  surface  temperature  T(o)  of  the  propellant  as  well  as  the  flux  q. 

The  flux  q  equals  qf  given  by  Equations  11  or  12  depending  upon  whether  or  not 
there  is  a  cross  flow.  The  surface  temperature  T(o)  may  be  calculated  in 
terms  of  time-dependent  fluxes  q  following  the  procedure  described  in  Appendix  A. 
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APPENDIX  E 


FLUID  DYNAMICS  AND  STRESS  WAVES 

This  appendix  presents  means  used  to  predict  gas  flows  wi thing  cracks, 
and  the  generation  and  effect  of  stress  wave  upon  cracks. 

1.  FLUID  DYNAMICS 

1.1  BASIC  EQUATIONS 

The  equations  of  fluid  dynamics  are  written  with  respect  to  the  Eulerian 
frame  of  reference  in  which  the  independent  variables  are  the  distance,  x, 
along  the  crack  and  the  time,  t.  These  equations  represent  the  conservation 
of  mass,  momentum,  and  energy.  The  present  gas  dynamic  model  includes  the 
effects  of  inertia,  wall  fraction,  mechanical  wall  response,  and  mass  and 
energy  addition.  Heat  transfer  effects  are  incorporated  into  the  energy 
addition  term. 

The  conservation  laws  can  be  developed  by  considering  a  control  surface 
enclosing  an  element  of  volume  of  length,  dx.  This  control  surface  is 
illustrated  in  Figure  E-l.  The  conservation  laws  state  that  the  time  rate  of 
change  of  the  entity  considered,  within  the  control  volume  must  equal  the  net 
flux  of  this  entity  through  the  control  surface,  plus  any  related  boundary 
contributions  such  as  mass  or  energy  addition,  work  done  at  the  boundary,  or 
boundary  forces.  The  area  of  the  flow  channel,  A,  is  one  of  dependent  variables 
to  be  considered.  However,  in  view  of  the  fact  that  the  lateral  extent  of 
a  crack  is  very  wide  compared  to  its  width  Cw.  We  shall  set  the  lateral  ex¬ 
tent  of  the  flow  channel  equal  to  unity  so  that  cross-section  area  A  of  the 
channel  is: 

A  -  1  •  Cw  (E-l) 

It  should  be  noted  that  although  this  variable  has  the  dimensions  of  length, 
a  dimensional  check  of  the  following  equations  will  imply  that  it  should  have 
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x  +  dx 


Figure  E- 


Control  surface  for  flow  channel 
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the  dimensions  of  length  squared.  This  is  the  case  when  it  represents  the 
area  of  the  flow  channel. 

The  three  conservation  equations  are: 


op  3p  •  d  u  3C  p  3C 

-_a  +  p  iiL  +  u  _£  =  JL  .  _ il  _£ _ w 

at  pg  3x  u  3x  cw  3*  cw  ^ 


(E-2) 


where: 


Pg  =  gas  density 
u  =  gas  velocity 

m  =  mass  addition  per  unit  length  per  unit  time 
The  wall  motion  (velocity)  is  given  by  ( 3Cw/3 t )  (see  Figure  E-l) 


Momentum 


9u  ,  3u  .  1  3P  2a  um 


(E-3) 


where  t  =  wall  shearing  stress.  The  wall  shearing  stress  can  be  expressed  in 
terms  of  the  conventional  pipe  friction  coefficient,  F,  defined  by: 


(E-4) 


Furthermore  this  stress  acts  in  a  direction  opposite  to  the  fluid  motion.  The 
momentum  equation  can  then  be  reformulated  as: 


+  u  +  JL  1?.  = 

at  u  ax  pg  ax 


(E-5) 


The  coefficient  F  is  assumed  to  be  constant  at  a  nominal  value  because  its 
value  and  dependence  are  uncertain  for  this  complex  flow  environment. 


EneAgy 


3e-r  3Cy  0 

IF  +  u  IF  +  FT  Tx  (Pu)  =  fV 


Q  Pu 

F "  F  '3x 


p  ’jw  v 
st  p9^ 


(E-'5) 
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where  Q  =  energy  addition  per  unit  Tenth  per  unit  time  and  ey  =  specific  total 
energy. 

The  specific  total  energy  is  given  by: 

eT  =  e  +  >su  (E-7) 

where  e  =  specific  internal  energy 

Equations  E-2,  E-5,  and  E-6  have  been  arranged  so  that  the  terms  which 
define  the  contributions  of  the  special  effects  are  grouped  on  the  right  hand 
side  of  the  equal  sign.  When  the  sum  of  the  terms  of  the  right  hand  side  are 
set  equal  to  zero  the  conventional  gas  dynamic  equations  for  the  nonsteady 
flow  in  a  constant  area  channel  are  obtained.  These  reduced  equations  will  be 
used  to  establish  the  numerical  method. 

1.2  EQUATION  OF  STATE 

The  equation  of  state  which  has  been  selected  initially  for  the  propellant 
reaction  products  is  that  of  a  perfect  gas.  This  equation  of  state  is  con¬ 
sidered  to  be  adequate  for  the  initial  phases  of  the  investigation,  but 
eventually  a  better  formulation  must  be  established. 

The  equation  of  the  state  of  the  gas  is  the  following: 

P  (1/Pg  -  b)  *  RT  (E-8) 

where  R  =  gas  constant  and  T  *  absolute  temperature  of  the  gas. 

Furthermore,  the  specific  heats  at  constant  pressure  and  at  constant 
volume  are  both  constant.  The  internal  energy  e  is  given  by: 

e  =  P  (1/Pg  -  b) / ( y  -  1)  (E-9) 

where  y  s  ratio  of  specific  heats. 

The  sound  velocity  of  the  gas,  c,  is  given  as: 

c  (E-10) 
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The  values  of  the  gas  constant,  the  two  specific  heats  and  the  ratio  of 
specific  heats  are  subject  to  some  uncertainty.  The  following  values  have 
been  selected  for  the  gas  constant  and  the  ratio  of  specific  heats  for  the 
reaction  products 

R  =  3228  (cm/°C) 

Y  =  1.2 

1.3  FLOW  CONFIGURATIONS  AND  BOUNDARY  CONDITIONS 

The  solutions  of  the  above  field  equations  are  subject  to  the  initial 

conditions  within  the  crack  and  to  boundary  conditions  at  both  the  upstream 

and  downstream  end  of  the  flow  channel.  The  initial  conditions  are  those  of 

a  gas  at  rest  (u(x)  =  0)  and  at  some  pressure  pQ,  an  temperature,  Tq.  The 

initial  height,  CwQ(x)  of  the  channel  is  also  specified,  however  no  Initial 

wall  motion  is  permitted  (i.e. ,  (ZCJ'dt)  =  0).  The  propellant  mass  in  the 

wu 

vicinity  of  the  crack  is  thus  in  mechanical  equilibrium  with  the  initial 
pressure  field  within  the  crack. 

Two  basic  crack  configurations  are  treated.  These  are  illustrated  in 
Figure  E-2.  They  consist  of  a  high  pressure  cavity  connected  to  the  upstream 
and  (x  =  0)  of  a  crack  of  length,  L  ,  together  with  one  of  two  downstream  end 
conditions.  One  configuration  has  a  simple  closed  end  while  the  other  con¬ 
figuration  consists  of  a  connection  to  a  low  pressure  cavity.  The  high 
pressure  cavity  contains  a  gas  at  a  pressure  pc  and  temperature  Tc,  both  of 
which  are  held  constant  with  respect  to  time.  The  gas  pressure  within  the 
low  pressure  cavity  is  held  constant  at  the  level  of  the  initial  gas  pressure 
within  the  crack  (i.e.,  PQ)  and  only  outflow  is  permitted. 

The  boundary  conditions  at  the  high  pressure  cavity  end  will  depend  upon 
whether  the  gas  flow  is  into  (inflow)  or  out  of  the  crack  (outflow).  These 
boundary  conditions  are  illustrated  in  the  Hodograph  plane  of  Figure  E-3.  The 
cavity  State,  Sc,  (a  rest  state)  is  the  appropriate  reference  state.  When 
the  pressure  is  low  relative  to  the  cavity  pressure  at  the  downstream  end  of 
the  crack  inflow  may  occur.  If  inflow  does  occur,  the  cavity  gas  will  flow 
into  the  crack  after  first  expanding  isentropical ly  (at  constant  energy)  In 
the  vicinity  of  the  inlet.  This  expansion  process  will  accelerate  the  gas 
and  lower  both  the  pressure  and  the  sound  velocity.  The  flow  Mach  number,  M 
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E"3,  Rodograph  plane  -  boundary  conditions  at  cavity  end. 


(M  =  u/c)  will  increase  until  it  reaches  a  state  which  is  in  equilibrium  with 
the  Internal  flow  at  the  boundary  provided  it  does  not  exceed  a  value  of 
unity.  At  that  critical  point  no  further  expansion  can  occur  since.  In  effect, 
no  further  Information  regarding  any  additional  expansion  can  be  communicated 
back  to  the  cavity.  Any  additional  expansion,  if  it  can  occur,  will  occur  as 
a  non-steady  expansion  within  the  crack.  The  field  equation  will  provide  for 
the  solution  of  this  additional  expansion.  Thus  the  permissible  boundary 
conditions  at  the  downstream  end  of  the  crack  under  Inflow  conditions  will  be 
those  states  associated  with  the  energy  line  bounded  by  the  cavity  State,  Sc> 
at  one  end  and  the  sonic  inflow  State  S$1.  at  the  other  end.  The  energy  line 
is  defined  by 


cc2  •  c2  +  u2  (E-ll) 

If  subsonic  outflow  occurs  the  pressure  at  the  upstream  end  of  the  crack 
will  be  equal  to  the  cavity  pressure,  Pc  (see  Figure  E-3).  However,  if  the 
flow  becomes  sonic  or  supersonic  then  the  pressure  can  change  to  any  value 
such  that  the  boundary  state  within  the  crack  lies  In  the  supersonic  outflow 
region  bounded  by  the  sonic  outflow  line.  This  region  is  illustrated  In 
Figure  E-3. 

The  same  boundary  conditions  apply  at  the  low  pressure  cavity  end  of  the 
crack  when  this  configuration  is  used,  however  the  present  model  does  not 
permit  inflow  to  occur.  Some  model  modifications  In  this  area  may  be 
needed  however,  they  will  be  influenced  by  the  physical  model  associated  with 
this  low  pressure  cavity,  such  as  for  example  a  cavity  filling  process  with 
a  subsequent  buildup  of  pressure  and  temperature  within  the  finite  volume 
cavity. 

The  boundary  condition  for  the  closed  end  is  a  simple  one  of  no  flow 
(i.e. ,  u(Lc)  =  0).  The  current  model  does  permit  the  crack  length  to  be 
extended  arbitrarily  whenever  this  configuration  is  used.  Whenever  the  crack 
Is  extended,  the  new  portion  Is  filled  with  a  rest  gas  at  the  initial  pressure 
and  temperature  (P0,  TQ)  and  the  closed  end  boundary  condition  is  applied  to 
the  new  closed  end  location.  This  crack  extension  is,  in  effect,  an 
Instantaneous  crack  extension. 
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1.4  NUMERICAL  PROCEDUIE 


The  numerical  solution  technique  employed  in  the  gas  dynamic  model  of  the 
present  study  is,  in  its  reduced  form,  *  rather  conventional  Eulerian  method 
of  the  FLIC  (fluid  in  cell)  type.  This  method  was  originally  developed  at  the 
Los  Alamos  Scientific  Laboratory. 

The  numerical  solution  of  the  foregoing  equations  proceeds  from  the 
initial  conditions  specified  in  a  forward  stepping  time  wise  manner  subject  to 
the  numerical  solution  to  the  field  equations,  auxilliary,  time  wise  Inputs, 
and  the  appropriate  boundary  conditions.  An  artificial  viscous  pressure 
term,  q,  is  added  to  the  thermodynamic  pressure  during  the  computations  when 
the  flow  is  subsonic  and  the  compression  rate  is  positive.  This  contributes 
to  the  suppression  of  flow  discontinuities  and  to  the  computational  stability 
in  regions  of  subsonic  flow. 

1.4.1  The  Computing  Me&k 

A  one-dimensional  mesh  of  uniform  length  ax,  cells  is  established.  Each 
cell  is  Identified  by  an  indice,  i,  corresponding  to  its  center.  Thus  the 
boundaries  of  the  i th  cell  are  at  the  location  i  t  Increasing  1  corres¬ 
ponds  to  increasing  x.  The  selection  of  uniform  cell  length  Is  an  Initial 
convenience.  If  subsequent  results  indicate  that  one  portion  of  the  crack 
requires  substantially  r.rrater  resolution  in  the  space  variable,  x,  than  do 
other  portions  then  variable  length  cells  can  be  introduced  with  little 
addi tonal  complications. 

The  values  of  the  dependent  variables  are  associated  with  the  cell  centers 
with  the  exception  of  the  pseudoviscous  terms  which  are  computed  at  the  cell 
boundaries.  Assuming  that  all  properties  are  known  for  each  cell  at  some  time 
tn,  the  computational  procedure  is  to  determine  tne  state  In  each  cell  at  a 
later  time  tn+l  =  tn  +  At.  The  time  step  At  is  restricted  in  magnitude  by 
conditions  required  for  stability  of  the  computations.  The  resulting  space 
time  grid  network  is  illustrated  in  Figure  E-4. 

1.4.2  Stabi&Lty  To/unuIntLon 

There  are  several  stability  type  conditions  that  over  Investigators  have 
found  applicable  with  this  type  numerical  technique.  One  such,  restriction  is 
that  (umavAt/Ax  <  1.  If  fluid  particles  were  explicitly  treated  In  the 

Ilia  a 
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computation,  this  criterion  would  prevent  a  particle  from  crossing  a  complete 
cell  during  one  time  step.  In  the  pure  Eulerian  scheme  the  Interpretation  is 
made  that  the  transport  terms  are  calculated  more  accurately  with  corresponding 
Improved  averaging  of  the  numerical  fluctuations.  A  second  restriction  Is  the 
familiar  Courant  condition,  cAt/Ax  <  1,  which  limits  the  propagation  of  fluc¬ 
tuations  in  the  subsonic  flow  regions.  Another  source  of  instability  Is  the 
computation  of  negative  Internal  energies  in  the  cell  due  primarily  to  the 
treatment  of  boundary  conditions,  but  also  from  Inherent  numerical  fluctuations, 
for  a  fluid  obeying  the  ideal  gas  relation,  the  criterion  that  prevents  nega¬ 
tive  Internal  energies  Is  that  4(y  -  1)  u„a„At/Ax  <  1.  These  stability 

lUOA 

criteria  can  be  written  as 


At  <  min  (£*■),  over  i  (E-12) 

ui 

At  <  min  (7^-) ,  over  i  ( E- 13 ) 

ci 

At  <  0.625  min  (7^-),  over  i  ( E- 14) 

ui 

A  simple  stability  criteria  which  satisfies  all  of  the  foregoing  requirements 
is 


4t  *  15  n'i"  (M»-(7r-,-c'TK  over  1  (E-15) 

where  the  factor  %  Is  chosen  to  reflect  the  typical  severity  with  which  the 
inequalities  (E-12)  to  ( E- 1 4 )  are  generally  applied.  More  details  regarding 
the  computations  procedure  are  described  In  Appendix  8  of  Reference  14. 

2.  GENERATION  AND  EFFECTS  OF  STRESS  WAVES 

The  flow  of  high  pressure  gases  Into  cracks  within  an  adjacent 
propellant  mass  and  any  subsequent  burning  of  these  propellant  surfaces  will 
create  a  time  varying  pressure  environment  with  the  cracks.  The  surrounding 
propellant  mass,  which  was  originally  in  mechanical  equilibrium  with  a  low 
pressure  distribution  within  the  crack  will  respond  to  these  new  mechanical 
loads.  It  Is  the  purpose  of  this  appendix  to  describe  several  mechanical 

IIT  RESEARCH  INSTITUTE 


95 


response  models  which  have  been  used  In  conjunction  with  the  gas  dynamic  model 
described  In  Section  1  of  this  appendix.  One  response  model  treats  the  early 
phases  of  the  mechanicl  response  when  stress  waves  radiate  from  the  crack  and 
the  crack  Is  in  a  growth  phase.  The  second  model  tests,  in  a  simple  fashion, 
the  Interaction  of  these  radiating  waves  with  the  confining  shell  structure 
and  their  subsequent  interaction  with  the  crack.  Under  these  conditions  the 
crack  growth  will  generally  be  arrested  and  may  ultimately  lead  to  partial  or 
complete  crack  collapse. 

2.1  MECHANICAL  PROPERTIES  OF  THE  PROPELLANT 

The  propellant  response  to  the  transient  pressure  loads  within  a  crack 
will  be  that  associated  with  the  generation  and  propagation  of  stress  waves 
in  the  propellant  mass.  The  mechanical  properties  of  most  propellant  materials 
are  quite  complex  in  that  they  are  generally  nonlinear,  rate  sensitive,  and 
hysteretic  in  nature.  Any  simple  mechanical  response  model  cannot  deal  effect¬ 
ively  with  these  types  of  complex  behavior  characteristics;  however  their 
impact  upon  the  subject  problem  is  not  considered  to  be  important.  Rather  it 
should  be  sufficient  to  model  in  some  simple  way  the  gross  compressibility 
characteristics  of  the  material.  Thus  the  material  will  be  viewed  Initially 
as  a  simple  linear  isotropic  elastic  material  which  does  not  yield  or  fail  in 
any  manner  under  the  imposed  stresses. 

An  elastic  material  will  support  two  basic  types  of  stress  waves,  i.e. , 
dilatation  and  shear  waves.  These  waves  will  propagate  at  the  following  speeds. 

V-clatcuUon  Wave.  Speed, 

cd  <£-16> 

SheaA  (itave  Speed,  c. 

tt  (E-17> 
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where: 


E'  =  Young's  modulus 
t  *  Poisson's  ratio 
p  «  densi  t.y 


The  ratio  of  these  two  wave  speeds  is  dependent  only  upon  Poisson's  ratio,  vis: 


c 

c 


£ 

d 


(E-18) 


The  influence  of  Poisson's  ratio  is  not  great  until  it  exceeds  a  value  of 
about  0.45  and  then  in  primarily  affects  the  dilitation  wave  speed.  The 
effective  wave  speed  c*  =  /t/p  is  the  most  significant  parameter.  Since 
Poisson's  ratio  may  be  in  the  broad  range  of  from  0.2  to  0.4  the  dilitation 
wave  speed  will  be  approximately  1.2  c*  while  the  shear  wave  speed  will  be 
approximately  0.6  c*.  Thus  the  wave  speed  ratio  will  be  approximately  0.5. 

The  density  of  the  propellants  are  well  known  and  a  nominal  value  of  1.5  g/cm3 
is  used  in  the  mechanical  response  models.  A  nominal  value  of  250  cm/ms  was 
selected  for  the  effective  wave  speed.  This  corresponds  to  a  value  of 
approximately  106  psi  for  Young's  modulus  and  yields  dilitation  and  shear 
wave  speeds  of  300  cm/ms  (approximately  10,000  fps)  and  150  cm/ms  respectively. 


2.2  INITIAL  RESPONSE  MODEL 

During  the  entry  phase  of  the  gas  flow  into  the  crack  a  shock  wave  is 
generated  which  propagates  at  a  speed  (depending  upon  a  number  of  parameters) 
of  approximately  150  cm/ms.  The  pressure  distribution  behind  this  shock  wave 
is  relatively  uniform.  Weak  disturbances  (i.e. ,  sound  waves)  propagate  at  the 
local  sound  speed  relative  to  the  gas.  The  sound  speed  is  generally  in  the 
approximate  range  of  from  50  to  100  cm/ms  while  the  particle  velocity  is  In 
the  range  of  from  -50  to  100  cm/ms.  Thus  these  disturbances  will  propagate, 
on  the  average,  at  an  absolute  velocity  of  about  100  cm/ms.  These  considera¬ 
tions  lead  us  to  the  following  approximate  inequality. 


|(u  i  c) |  <  U  «  c$  <  cd 


(E-19) 
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where: 


U  =  shock  velocity  in  the  gas 
m  =  particle  velocity  in  the  gas 
c  =  sound  velocity  in  the  gas 

Thus  the  wave  system  illustrated  in  Figure  E-5  should  exist  initially. 
While  this  wave  system  is  quite  complex,  many  of  the  waves  should  be  weak  and 
need  not  be  considered.  This  is  especially  true  for  the  outrunning  waves. 

The  primary  motion  of  the  propellant  mass  will  be  in  a  direction  normal  to  the 
crack,  causing  the  crack  width  to  increase  locally.  For  this  reason,  a 
very  simple,  one  dimension  wave  propagation  (plane  strain)  model  has  been 
selected  with  which  to  define  the  response  of  the  propellant  adjacent  to  the 
crack.  The  momentum  equation  for  this  case  is  of  the  form 

At  =  (ogcd)  (E-20) 


where: 

At  =  change  in  stress  at  the  crack  boundary 
Au  =  change  in  velocity  at  the  crack  boundary 

The  change  in  the  stress  at  a  given  location  along  the  crack  is  identical  to 
the  change  in  the  corresponding  gas  pressure.  Thus  the  wall  velocity  W  used 
in  the  gas  dynamic  model  is  given  for  the  ith  cell  as 


wi  3  <<*«,)  <Pi  '  Po) 


(E-21) 


where: 

pcd  =  shock  impedence  (=450  bars-ms/cm) 

P.j  =  gas  pressure 

PQ  =  Initial  gas  pressure 

This  expression  was  modified  to  include  some  contribution  from  the  adjacent 
cells  (1-1,  and  i  +  1).  In  this  manner  some  influence  of  the  spreading  nature 
of  the  propagation  within  the  propellant  mass  was  included,  at  least  in  a 
crude  manner. 
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2.3  INFLUENCE  OF  REMOTE  BOUNDARIES 


The  stress  wave  which  radiates  away  from  the  crack  will  eventually  react 
with  the  remote  boundaries  of  the  propellant  mass,  reflect  and  return,  in  some 
modified  form,  to  the  crack.  These  reflected  waves  will  also  influence  the 
response  crack  walls.  This  type  of  wave  system  is  illustrated  in  part  a  of 
Figure  E-6.  The  nature  of  reflected  wave  system  will  be  very  complex  and 
varied.  For  this  reason  a  number  of  idealized  reflection  models  will  have  to 
be  established  to  define  the  configurations  and  conditions  of  interest. 

As  an  initial  attempt  to  treat  the  influence  of  the  remote  boundaries  of 
the  propellant  mass  one  idealized  model  was  established.  This  model  is 
shown  in  part  b  of  Figure  E-6.  It  is  designed  to  introduce  the  ultimate 
confining  influence  of  a  cased  mass  of  propellant.  A  radiative  stress  field, 
APr,  which  is  the  integrated  or  averaged  value  of  the  current  local  stress 
field  along  the  entire  crack  is  defined  and  assumed  to  propagate,  under  the 
conditions  of  the  plane  strain,  into  the  propellant  mass.  These  stress  waves 
interact  with  the  remote  boundary  after  a  delay  time  corresponding  to  a 
boundary/crack  separation  distance,  Lg.  The  case  is  treated  as  a  spring 
supported  mass  characterized  by  its  inertia  (weight  per  unit  area)  and 
stiffness  (breathing  mode  of  the  case).  The  reflected  wave  system,  AP$  is 
evaluated  from  this  case  boundary  interaction  and  applied  uniformly  along  the 
crack  after  the  appropriate  delay  period.  This  reflected  wave  system  then 
interacts  locally  with  the  crack  boundary  such  that  the  wall  velocity  changes 
according  to  the  following  momentum  consideration 

AWi  =  -2APs/(pC$)  (E-22) 
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E-6.  Influence  of  remote  boundaries  on  mechanical-  response. 


APPENDIX  F 

NOMENCLATURE  AND  HMX  PROPERTY  DATA 


This  appendix  described  variables  used  in  this  report  as  well  as  the 
values  used  for  the  properties  of  HMX  and  Its  reaction  gases. 

1.  NOMENCLATURE 

Nomenclature  Is  presented  below.  In  the  report,  bars  over  the  parameters 
indicate  steady-state  values  while  the  subscript  o  indicates  initial  values  or 
values  ironediately  prior  to  the  event  under  consideration. 

a  rate  of  steady  burning  at  ambient  pressure,  cm/sec. 

A  area  of  flow  channel,  cm  (see  Appendix  E  for  description  of  units). 

A^  bearing  area  between  piston  and  chamber,  cm2. 

A  area  of  piston  face,  cm2, 

p 

B  constant  used  to  evaluate  erosive  heating. 

c  sound  speed,  cm/msec. 

cd  dill tati on  wave  speed,  cm/msec. 

c$  shear  wave  speed,  cm/msec. 

C  constant  given  by  2  J/KpC  ,  cm2  °K  (sec)0 ' 5/cal . 

C  specific  heat  of  gases  evolved  by  propellant  at  constant  pres- 

9  sure,  cal/g-°K. 

Cm  specific  heat  of  moltens  propellant,  cal/g-°K. 

Cp  specific  heat  of  solid  propellant,  cal/g-°K. 

C  crack  width,  cm. 

W 

ci,  cz  constants  used  to  describe  heat-transfer  coefficient  h  where  h  = 
Cl[Z  exp  <-E/Tf)]C2 

d  drag  coefficient 

Dp  piston  travel  needed  to  close  vent  holes,  cm. 

i  internal  energy  of  combustion  gases,  cal/g. 

E  activation  energy  of  propellant  dividied  by  gas  constant,  °K. 

E1  Young's  modulus,  psi. 

f  rates  of  crack  area  to  that  of  planar  surface,  dimensionless. 

F  conventional  pipe  friction  coefficient, 

g  acceleration  of  gravity. 
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heat-transfer  coefficient,  cal/cm2-sec-°K. 
convective  heat-transfer  coefficient,  cal/cm2-sec-°K. 
subscript  indicating  value. 

mechanical  impedence  of  solid  propellant,  equals  I  +  0.002P, 
bars-sec/cm. 

see  I. 

subscript  Indicating  value  of  parameter  during  time  step  t^. 
thermal  conductivity  of  solid  propellant,  cal/cm-sec-°K. 

Indepth  crack  distance,  cm. 
mass,  g. 

mass  of  driver,  lbs. 

mass  of  unit  area  of  foam  layer,  g/cm2. 

mass  of  combustion  gases. 

mass  of  piston  assembly,  lbs. 

molecular  weight  of  reaction  gases,  g/mole. 

exponent  of  pressure  P  used  to  describe  steady  burning  rate  r^. 

pressure,  bars  or  k  bars. 

maximum  value  of  P  achieved,  k  bars. 

Prandth  number,  dimensionless. 

anplitude  of  incident  stress  wave,  bars. 

heat  flux  entering  solid  propellant,  cal/cm2-sec. 

rate  of  heating  of  unit  area  of  foam,  cal /cm2 -sec. 

heat  flux  applied  at  depth  x*j  during  time  step  Atj,  cal/cm2  -sec. 

heat  flux  from  foam  to  melt  interface,  cal/cm2-sec. 

mean  conductive  heat  flux  q^,  cal/cm2-sec. 

sensible  heat  per  unit  mass  of  foam,  cal/g. 

latent  heat  of  fusion  of  propellant,  cal/g. 

reaction  heat  of  propellant,  cal/g. 

heat  generated  within  foam  per  unit  mass  of  evolved  propellant 
gas,  cal/g. 

rate  of  melting  of  propellant,  cm/sec. 
rate  of  propellant  burning,  cm/sec. 
gas  constant,  cm/°K. 
time,  sec. 

j  th  time  step,  sec. 
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gas  velocity,  cm/msec. 

shock  velocity,  cm/msec. 

temperature  of  foam  or  propellant  surface,  °K. 

temperature  of  combustion  gas,  °K. 

melt  temperature  of  propellant,  °K. 

Initial  temperature  of  propellant,  °K. 
velocity  of  driver,  cm/sec. 
velocity  of  piston,  cm/sec. 

velocity  of  crack  wall  exposed  to  Incident  stress  wave,  cm/sec. 

velocity  of  crash  wall  not  exposed  to  Incident  stress  wave, 
cm/sec. 

distance,  cm. 

depth  of  melt  Interface  at  end  of  time  step  Atj,  cm. 

depth  at  which  flux  q1^  is  applied  during  time  step  At^,  cm. 

frequence  factor  associated  with  propellant,  1/sec. 

thermal  dlffusivity  of  propellant,  cm2/sec. 

dimensionless  term  used  to  determine  time  steps  used  to  calcu¬ 
late  propellant  temperature. 

ratio  of  specific  heats  of  combustion  gases,  dimensionless, 
particle  velocity  of  gas,  cm/msec. 

factor  used  to  determine  depths  x'.  at  which  fluxes  q1.  are 
applied,  dimensionless.  J  J 

density  of  solid  propellant,  g/cm* . 

density  of  combustion  gases,  g/cm* . 

Poisson's  ratio. 


2.  HMX  PROPERTY  OATA 

Properties  for  HMX  propellant  and  evolved  gases  are  presented  In 
Table  F-l. 
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TABLE  F-l.  PROPERTIES  OF  HHX  PROPELLANT  AND  GASES 


Parameter 

Value 

Source 

Constants  a,  n  of  Equation  (2) 

a  *  0.030  cm/sec 
n  =  0.86  (dimensionless) 

Reference  7 
Reference  7 

Constants  Cj ,  c2  of  Equation  (9) 

ci  =  7.6*10®  cal/cm2-sec-°K 

See  Appendix  B 

C9 

c2  *  0.338  (dimensionless) 
0.5  cal/g-°K 

See  Appendix  B 
Assumed 

Cp 

0.4  cal/g-°K 

Reference  1 

Sn 

0.5  ca1/g-°K 

Reference  1 

f 

1.5  (dimensionless) 

Assumed 

E 

27, 000° K 

Reference  13 

K 

0.0013  cal/cm-sec-°K 

Assumed 

% 

50  cal/g 

Reference  1 

<V 

3517  cm/°K 

Assumed 

R 

3517  cm/°K 

Assumed 

Tm 

555°  K 

Reference  1 

2 

0.5  10*%/sec 

Reference  13 

p 

1.9  g/cm 

Reference  7 

Y 

1.2 

Assumed 
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